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Abstract 


Myers,  Clifford  A.,  Carleton  B.  Edminster,  and  Frank  G.Hawksworth. 
1976.  SWYLD2:  Yield  tables  for  even-aged  and  two-storied  stands 
of  southwestern  ponderosa  pine,  including  effects  of  dwarf 
mistletoe.  USDA  For.  Serv.  Res.  Pap.  RM-163,  25  p.  Rocky 
Mt.  For.  and  Range  Exp.  Stn.,  Fort  Collins,  Colo.  80521. 

Presents  a  procedure  for  computation  of  yield  tables  for  pon- 
derosa pines  in  Arizona  and  New  Mexico.  Possible  alternatives 
include:  even-aged  or  two-storied  stands,  healthy  or  diseased  stands, 
and  managed  or  unmanaged  stand  densities.  Stand  conditions  and 
severity  of  dwarf  mistletoe  infestation  change  with  time  and  in 
response  to  intermediate  cuttings.  Supersedes  SWYLD,  published 
in  1972  as  USDA  For.  Serv.  Res.  Pap.  RM-87.  A  concise  user's  guide 
for  program  SWYLD2  is  available  as  USDA  For.  Serv.  Gen.  Tech. 
Rep.  RM-23,  1976. 

Keywords:  Stand  yield  tables,  timber  management,  forest  manage- 
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Introduction 

Computer  program  SWYLD  (Myers  and  others 
1972)  was  developed  to  compute  yield  tables  for 
managed,  even-aged  stands  of  southwestern  pon- 
derosa pine  (Pinus  ponderosa  Laws.).  Users  of 
SWYLD  recognized  a  need  for  a  capability  of 
handling  a  wider  range  of  situations.  A  few  of  the 
original  relationships  and  analytic  procedures 
have  been  modified  and  new  ones  have  been  added. 
Program  SWYLD2,  described  in  this  paper,  has 
been  written  to  supersede  SWYLD,  published  as 
a  USDA  Forest  Service  Research  Paper  RM-87. 

Modifications  and  additions  now  provide  the 
following  features:  (1)  yield  tables  are  computed 
for  even-aged  stands  or  for  each  story  of  two- 
storied  stands;  (2)  equations  predict  the  contribu- 
tion of  the  overstory  to  dwarf  mistletoe  infection 
in  the  understory,  if  any;  (3)  height  and  diameter 
growth  equations  apply  to  a  wide  range  of  stand 
densities,  and  (4)  an  intermediate  or  regeneration 
cut  is  no  longer  required  at  the  initial  age  printed 
in  the  yield  table. 

The  material  presented  includes:  (1)  descrip- 
tions of  the  field  measurements,  (2)  descriptions 
of  the  relationships  obtained  from  the  measure- 
ments, (3)  a  computer  program  that  contains  the 
relationships,  and  (4)  an  example  of  what  the  pro- 
gram can  produce.  Instructions  for  adapting  the 
procedures  and  program  to  other  species  or  forest 
regions  are  also  given. 

A  published  summary  of  instructions  on  prep- 
aration of  yield  tables  for  diseased  stands,  di- 
rected to  personnel  of  pest  detection  and  control 
programs,  is  mostly  applicable  to  SWYLD2 
(Hawksworth  and  Myers  1973).  One  important 
exception  is  that  the  new  program  is  useful  for 
two-storied  stands  as  well  as  those  that  are  even- 
aged.  A  second  exception  is  that  SWYLD2  accepts 
more  than  one  way  of  expressing  degree  of  dwarf 
mistletoe  infestation. 


SWYLD2  is  written  in  standard  FORTRAN 
IV  and  can  be  run  on  almost  any  computer  that 
provides  50,000  octal  memory  locations.  Yield 
tables  can  be  produced  at  a  cost  of  about  8  cents 
each,  excluding  cost  of  program  compilation. 

Uses  of  Yield  Tables 

Yield  tables  are  essential  guides  for  forest 
managers.  They  report  probable  wood  yields  that 
result  from  specified  combinations  of  such  factors 
as  site  quality,  utilization  standards,  and  fre- 
quency and  intensity  of  thinning.  They  are,  there- 
fore, the  basis  for  timber  management  planning. 
They  also  provide  an  important  part  of  the  infor- 
mation needed  for  determining  the  influence  of 
timber  treatments  on  other  forest  resources. 

Two  types  of  yield  tables  can  be  computed  by 
program  SWYLD2:  (1)  tables  for  managed  stands, 
and  (2)  tables  that  describe  the  remainder  of  the 
rotation  or  life  of  existing  stands.  The  term  "man- 
aged, ' '  as  used  in  this  paper,  means  stand  density 
is  controlled  throughout  the  life  of  the  stand.  The 
term  "stand,"  following  the  usual  definition, 
refers  to  an  area  of  trees  consistent  in  site  quality, 
species,  structure,  density,  and  incidence  of  dis- 
ease. Uniformity  of  site  quality  and  stand  char- 
acteristics must  be  such  that  site  index  and  other 
area  averages  really  describe  the  unit  of  forest 
simulated.  Stand  area,  therefore,  may  range  from 
a  very  few  acres  in  rough  terrain  to  a  hundred  or 
more  acres  where  meaningful  averages  may  be 
computed. 

Managed  Stands 

Yield  tables  for  managed  stands  are  useful 
regardless  of  the  current  level  of  management. 
Well-managed  forests  can  benefit  from  refine- 


merits  in  operations  that  are  guided  by  compari- 
sons of  actual  conditions  with  a  good  standard. 
Where  conversion  to  managed  stands  is  under- 
way, yield  tables  provide  goals  toward  which 
conversion  can  be  directed. 

Several  factors  influence  the  ways  that  yield 
tables  should  be  produced  and  the  types  of  data 
used  to  compute  them.  These  factors  include: 

1.  Yield  tables  for  managed  stands  are  needed 
even  for  regions  and  species  where  managed 
stands  do  not  yet  exist.  Temporary  plots  must  be 
used  to  obtain  needed  relationships. 

2.  No  single  yield  table  can  be  adopted  as 
standard  for  a  species  or  region.  Applied  forest 
management  requires  the  use  of  various  practices 
that  can  and  will  differ  from  one  forest  to  another, 
and  must  appear  as  variables  in  yield  table  com- 
putations. 

3.  A  manager  should  not  be  restricted  to  only 
one  yield  table  per  working  group,  or  series  of 
stands  managed  under  the  same  silvicultural  sys- 
tem. He  must  have  the  opportunity  to  examine 
the  probable  future  results  of  current  operations, 
to  propose  changes  in  the  management  of  re- 
sources in  his  care,  and  to  estimate  the  effects  of 
these  changes  before  money  is  spent  on  them. 

The  second  and  third  items  above  relate  to  a 
disadvantage  of  normal  and  empirical  yield  tables. 
One  table  can  report  probable  yields  for  only  one 
combination  of  stand  characteristics  and  manage- 
ment objectives.  This  is  not  an  effective  basis  for 
comparing  alternatives  and  making  decisions. 
Since  a  group  of  yield  tables  can  show  various 
combinations  of  characteristics  and  objectives,  a 
means  of  producing  many  yield  tables  at  low  cost 
can  be  useful  to  a  forest  manager.  A  computer 
program  that  will  do  this  is  described  below  and 
listed  in  appendix  1.  Once  the  relationships  be- 
tween stand  variables  have  been  established,  a 
manager  can  examine  the  probable  outcomes  of 
many  possible  variations  in  management.  There 
is  no  need  to  delay  decisions  or  to  speculate  on 
what  may  happen  if  a  condition  or  procedure 
changes. 

Existing  Stands 

Existing  stands  are  mentioned  separately 
from  managed  stands  to  emphasize  that  yield 
tables  have  uses  besides  those  for  setting  man- 
agement goals.  In  this  case,  the  question  is  not 
the  proper  management  for  a  specific  combination 
of  tree  species,  region,  and  objective.  Rather,  the 
problem  is  to  select  the  proper  procedure  for 
changing  a  stand  condition  that  actually  exists 


to  some  desired  condition.  Change  may  be  called 
for  because  stand  density  is  too  high,  tree  ages 
are  greater  than  planned  felling  age,  or  disease 
problems  prevent  application  of  usual  manage- 
ment practices. 

As  an  example  of  the  problems  created  by 
existing  unsatisfactory  stand  conditions,  consider 
a  parasite  or  other  disease  on  ponderosa  pines. 
Management  decisions  must  be  made  for  each 
diseased  stand.  Possible  alternative  treatments 
are:  (1)  sanitation  thinning,  (2)  removal  of  a  dis- 
eased overstory  to  release  a  well-stocked  under- 
story,  (3)  no  treatment,  with  the  hope  that  the 
stand  may  eventually  become  merchantable,  or 
(4)  destruction  and  regeneration.  Knowledge  of 
stand  structure  and  density  will  remove  some 
alternatives  from  further  consideration.  For  fur- 
ther selection,  an  estimate  of  potential  yield  with 
each  alternative  is  necessary.  This  information 
can  be  provided  by  yield  tables  that  describe  not 
managed  or  goal  conditions,  but  possible  results 
from  treating  a  situation  that  actually  exists.  The 
final  alternative  chosen  must  then  also  consider 
other  forest  resources,  including  scenic  values. 

Importance  of  Dwarf  Mistletoe 

Yield  tables  for  healthy  stands  do  not  provide 
all  the  information  needed  for  decisionmaking  in 
the  Southwest.  For  practical  application,  impor- 
tant and  predictable  causes  of  reduced  growth 
and  yield  must  be  included  in  the  computations. 
Dwarf  mistletoe  {Arceuthobium  vaginatum  subsp. 
cryptopodum  (Engelm.)  Hawksw.  &  Wiens)  is  so 
widespread  and  so  reduces  growth  and  yield  that 
it  is  an  important  factor  in  management  planning. 
Andrews  and  Daniels  (1960)  found  this  parasite 
on  36  percent  of  some  2,700  ponderosa  pine  plots 
located  throughout  Arizona  and  New  Mexico.  On 
the  Lincoln  National  Forest  and  adjacent  Mes- 
calero  Apache  Reservation  in  southern  New 
Mexico,  more  than  50  percent  of  the  plots  were 
infested  (Hawksworth  and  Lusher  1956,  Andrews 
and  Daniels  1960). 

Dwarf  mistletoe  is  one  of  the  four  major 
causes  of  mortality  in  southwestern  ponderosa 
pine  (Pearson  1939,  Myers  and  Martin  1963). 
Pearson  (1938)  found  mortality  in  heavily  infested, 
cutover  stands  to  be  about  five  times  that  in 
comparable  lightly  infested  or  healthy  stands.  On 
the  Mescalero  Apache  Reservation,  mortality  in 
stands  with  dwarf  mistletoe  was  nearly  twice  as 
high  as  in  healthy  stands  (Hawksworth  and 
Lusher  1956).  The  differences  were  most  pro- 
nounced in  cutover  areas,  where  the  mortality 
rate  in  infested  stands  was  3.3  times  that  in 
mistletoe-free  stands. 


Several  studies  have  documented  the  adverse 
effects  of  dwarf  mistletoe  on  the  height  and  diam- 
eter growth  of  individual  ponderosa  pines  (Korstian 
and  Long  1922,  Sperry  1934,  Hawksworth  1961). 
Such  studies  do  not,  however,  reveal  the  total 
stand  loss  due  to  dwarf  mistletoe.  They  do  not 
report  mortality,  but  give  information  on  surviv- 
ing trees  only. 

A  recent  study  compared  growth  rates  and 
mortality  in  treated  and  untreated  mistletoe- 
infested  stands  in  Arizona  (Lightle  and  Hawks- 
worth  1973).  Some  of  the  information  in  SWYLD2 
on  the  effects  of  dwarf  mistletoe  was  obtained 
from  the  plots  used  for  that  study. 

The  publications  cited  above  dealt  primarily 
with  the  cumulative  effects  of  dwarf  mistletoe  in 
unmanaged  stands.  They  do  not  provide  data  that 
can  be  used  directly  for  yield  prediction  in  man- 
aged stands  or  for  comparisons  of  alternatives. 
Many  of  the  pine-mistletoe  relationships  needed 
for  these  purposes  were  reported  in  the  publica- 
tion describing  program  SWYLD  (Myers  and 
others  1972).  Modifications  of  some  of  these  and 
additional  relationships  are  reported  here  for  the 
first  time. 

The  measure  of  dwarf  mistletoe  infestation 
used  in  SWYLD2  is  the  6-class  rating  system 
developed  by  Hawksworth  (1961).  A  rating  is 
obtained  for  each  tree  (fig.  1).  The  individual  rat- 
ings are  then  averaged  to  obtain  the  value  for  the 
entire  stand  that  is  used  in  the  relationships  of 
program  SWYLD2. 


Description  of  Program  SWYLD2 

Program  SWYLD2  consists  of  a  main  pro- 
gram and  eight  subroutines.  Operations  performed 
by  each  routine  are  described  below  and  identified 
by  comment  statements  in  the  listing  of  the  source 
program  (appendix  1).  Subsequent  sections  of  this 
paper  describe  the  relationships  used  to  compute 
values  of  the  dependent  variables.  Input  variables 
and  data  deck  structure  are  described  in  the  sec- 
tion headed  User-Supplied  Information. 


Main  Program 

The  main  program  performs  the  following 
seven  operations: 

1.  Reads  two  data  cards  containing  the  num- 
ber of  stands  to  be  processed,  the  length  of  the 
projection  period  in  years  (10.0  for  SWYLD2), 
and  minimum  commercial  cuts  in  board  and  cubic 
feet. 

2.  Checks  that  the  number  of  stands  to  be 
processed  is  not  zero  or  blank.  An  error  message 
is  printed  and  the  entire  job  terminates  if  this  error 
is  found. 

3.  Calls  subroutine  BEGIN  at  the  beginning 
of  processing  of  each  stand  to  set  values  of  many 
variables  to  zero  and  to  read  all  data  applicable 
to  the  stand  being  processed. 


INSTRUCTIONS 

STEP  1.  Divide  live  crown  into  thirds 

STEP  2.  Rate  each  third  separately. 
Each  third  should  be  given  a  — 
rating  of  0,  1  or  2  as  described 
be  1 ow. 

(0)  No  visible  infections. 

(1)  Light  infection  (1/2  or 
less  of  total  number  of 
branches  in  the  third  infected) 

(2)  Heavy  infection  (more 
than  1/2  of  total 
number  of  branches  in 
the  th  i  rd  i  nfected) . 

STEP  3.  Final ly,  add 
ratings  of  thi  rds  to 
obtain  rating  for 
total  tree. 


EXAMPLE 

If  this  third  has  no  visible 
infections,   its  rating  is  (0) 


f  this  third  is  lightly  infected, 
its  rating  is  (1 ) . 


f  this  thi  rd  is  heavi ly 
infected,   its  rating  is  (2). 


The  tree  in  this  example 
will  receive  a  rating  of 
0+1+2=3- 


Figure  1.— Instructions  for  and  example  of  the  use  of  the  6-class  mistletoe  rating  system  (Hawksworth  1961). 
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4.  Checks  many  variables  of  the  data  deck  for 
unwanted  zero  or  blank  values.  An  error  message 
is  printed  and  processing  terminates  for  any  stand 
with  such  errors  in  the  data  cards.  Understory 
variables  are  not  checked  because  they  must  have 
zero  or  blank  values  for  single-storied  stands. 

5.  Calls  subroutine  FIRST  to  compute  initial 
values  of  several  stand  characteristics  and  pro- 
gram controls. 

6.  Calls  2  subroutines  (CUTS,  PROJ)  in 
sequence  as  many  times  as  necessary  to  simulate 
the  number  of  projection  periods  required  to  reach 
the  final  age  desired  for  each  yield  table. 

7.  Calls  subroutine  TABLE  after  each  stand 
is  processed  to  print  yield  tables  for  the  stand. 
One  table  is  produced  for  single-storied  stands. 
Separate  tables  are  printed  for  the  overstory  and 
understory  of  two-storied  stands. 


Subroutine  BEGIN 

Subroutine  BEGIN  is  called  by  the  main  pro- 
gram to  initiate  the  processing  of  each  set  of  stand 
data.  Many  variables  are  assigned  initial  values 
of  zero  to  clear  computer  memory  of  data  from 
previous  computations,  and  to  prevent  processing 
of  an  understory  in  single-storied  stands.  Four 
cards  (types  3,  4,  5,  and  6)  are  read  by  BEGIN  to 
describe  stand  conditions  and  the  management 
controls  to  be  used  in  computations.  These  stand 
and  control  variables  are  described  in  a  subsequent 
section  and  in  the  listing  of  the  Order  and  Con- 
tents of  the  Data  Deck. 

Values  of  21  input  variables  are  stored  in  an 
array  named  VAR.  Some  items  in  VAR  are  ex- 
amined by  the  main  program  to  locate  errors,''  as 
already  described.  Others  are  used  by  other  sub- 
routines as  a  record  of  initial  values. 

The  data  deck  for  one  job  must  contain  one 
set  of  four  cards  (types  3, 4,  5,  and  6)  for  each  stand 
to  be  processed.  The  number  of  sets  or  stands  is 
specified  by  a  value  read  by  the  main  program. 


Subroutine  FIRST 

FIRST  computes  values  of  several  variables 
that,  with  the  input  variables,  describe  the  stand 
and  the  management  decisions  applicable  to  it. 
Where  appropriate,  separate  values  are  computed 
for  each  story  of  two-storied  stands.  The  final  age 
to  appear  in  the  yield  table  and  the  initial  basal 
area  are  computed.  Initial  dwarf  mistletoe  rating 
and  average  heights  of  dominant  and  codominant 
trees  are  also  computed  if  these  values  have  not 
been  entered  for  the  stand. 


Subroutine  SWVOL  is  called  by  FIRST  to 
compute  initial  stand  volume  per  acre  in  each  of 
the  units  provided  for  in  SWVOL.  These  volumes 
plus  other  computed  values  are  then  stored  for 
later  recall  when  the  yield  tables  are  printed. 


Subroutine  CUTS 

CUTS  executes  the  thinnings  and  regenera- 
tion cuts  called  for  by  the  variables  specifying 
management  controls.  CUTS  also  makes  any 
necessary  changes  in  the  interval  between  cuts 
and  in  residual  density  of  shelterwoods. 

A  scheduled  thinning  will  be  bypassed  if  any 
one  of  the  three  following  conditions  exists  at  the 
time  the  thinning  should  occur: 

1.  Comparison  of  current  growing  stock  level 
with  the  growing  stock  goal  indicates  that  the 
stand  or  story  is  already  understocked  in  relation 
to  the  goal. 

2.  Current  dwarf  mistletoe  rating  exceeds  the 
maximum  value  specified  by  the  program  user. 
This  value  will  usually  be  3.0  for  the  host-parasite 
combination  of  SWYLD2. 

3.  Noncommercial  volumes  would  be  cut  and 
the  appropriate  control  variable  is  given  a  value 
to  deny  this  option. 

Each  partial  cut  executed  is  followed  by  com- 
putation of  new  values  for  each  average  stand 
d.b.h.,  average  stand  height,  and  dwarf  mistletoe 
rating.  Each  new  diameter  and  height  is  larger  or 
smaller  than  the  corresponding  old  value,  depend- 
ing on  the  type  of  cutting  executed.  Post-cutting 
basal  areas,  residual  volumes,  and  volumes  re- ; 
moved  are  also  computed  and  stored  for  later 
printing  in  the  yield  tables. 

In  two-storied  stands,  CUTS  executes  thin- 
nings and  regeneration  cuts  separately  for  each 
story.  The  user  may  specify  simultaneous  cuttings 
in  both  the  overstory  and  understory  or  a  cutting 
in  only  one  story.  Growing  stock  level  specified 
for  a  cut  applies  only  to  the  basal  area  and  average 
diameter  of  the  story  involved  and  not  to  the  com- 
bined basal  area  of  both  stories. 

The  first  thinning  in  diseased  stands  will  be 
from  above,  as  simulated  by  subroutine  SWCUT2, 
unless  the  option  of  preventing  thinning  from 
above  is  exercised.  Average  diameters  and  heights 
will  therefore  be  reduced.  Subsequent  thinnings 
in  diseased  stands  will  increase  average  d.b.h.  and 
height,  but  by  lesser  amounts  than  in  healthy 
stands,  where  the  smaller  trees  make  up  a  larger 
percentage  of  those  removed.  This  effect,  an  alter- 
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native  computation  in  SWCUT1,  has  been  ob- 
served in  subsequent  thinnings  of  actual  diseased 
stands. 


Subroutine  PROJ 

Subroutine  PROJ  contains  the  relationships 
of  SWYLD2  that  describe  growth  and  other 
changes  over  time.  The  routine  consists  of  one 
long  loop  to  process  the  values  of  each  story  of 
two-storied  stands  separately. 

Values  computed  by  PROJ  during  each  pro- 
jection period,  if  needed,  are  as  follows:  (1)  non- 
zero dwarf  mistletoe  rating  (DMR)  that  accounts 
for  latent  infestation  after  apparent  elimination 
of  the  parasite  by  partial  cutting,  (2)  initial  DMR 
if  infestation  occurs  after  the  simulation  period 
begins,  (3)  value  of  DMR  at  the  end  of  the  period 
if  an  infestation  exists  and  if  neither  of  the  pre- 
vious situations  apply,  (4)  actual  growing  stock 
level,  (5)  average  d.b.h.  and  height  at  end  of  the 
period,  and  (6)  mortality  during  the  period.  Com- 
putations of  d.b.h.,  height,  and  mortality  include 
the  effects  of  dwarf  mistletoe,  if  present.  PROJ 
I  uses  the  combined  basal  area  of  both  stories 
(SBAS)  in  calculations  of  growing  stock  level  for 
each  story,  growth  projections,  and  mortality 
calculations. 

PROJ  completes  the  operations  of  each  pro- 
jection period  by  calling  SWVOL  to  compute 
volumes.  These  and  other  values  are  then  stored 
for  printing  at  the  end  of  the  simulation. 

Subroutine  TABLE 

TABLE  is  called  by  the  main  program  after 
each  set  of  stand  data  is  processed.  In  two-storied 
stands,  yield  tables  for  the  overstory  and  under- 
story  are  printed  separately  and  labeled  to  identify 
the  stand  component  described.  Only  one  table, 
without  a  story  designation,  is  printed  for  single- 
storied  stands.  Examples  of  the  tables  are  given 
in  appendix  2. 

Footnotes  to  each  table  report:  (1)  merchant- 
ability limits  for  cubic-  and  board-foot  volumes, 
(2)  type  of  initial  thinning  allowed,  (3)  minimum 
DMR  that  prevents  intermediate  cuts,  (4)  mini- 
mum volumes  included  in  total  yields,  (5)  controls 
on  noncommercial  thinnings,  (6)  initial  status  of 
dwarf  mistletoe  infestation  if  DMR  is  not  read 
from  data  card  type  5,  and  (7)  skipping  of  sched- 
uled intermediate  cuts. 


Subroutine  SWVOL 

Subroutine  SWVOL  is  called  by  other  sub- 
routines after  each  growth  projection  and  after 
each  intermediate  or  regeneration  cut.  Total  cubic 
volume  per  acre  (TCF)  of  all  trees,  from  groundline 
to  tip,  is  computed  with  stand  volume  equations. 
Factors  are  then  computed  and  used  to  convert 
total  cubic  volumes  to  merchantable  cubic  feet 
(FCTR)  and  to  board  feet  (PROD).  Utilization 
standards  for  these  units  are  given  in  the  comment 
statements  of  SWVOL  and  in  footnotes  to  the 
yield  tables. 

Two  or  three  equations  each  are  used  to  com- 
pute TCF,  FCTR,  and  PROD  to  avoid  the  need 
for  complex  forms.  These  relationships  are  not 
linear  over  the  necessary  ranges  of  the  indepen- 
dent variables. 

Conversions  to  other  units  or  utilization  stan- 
dards may  supplement  or  replace  those  already 
in  SWVOL  (appendix  1).  Additional  units  might 
be  square  feet  of  veneer  or  wood  weight  in  pounds 
(Myers  1960). 

Subroutine  SWCUT1 

Subroutine  SWCUT1  computes  average  stand 
d.b.h.  after  partial  cuttings  that  emphasize,  but 
are  not  restricted  to,  removal  of  the  smaller  trees. 
Average  diameter  and  height  after  cutting,  there- 
fore, are  larger  than  corresponding  values  before 
cutting.  The  percentage  of  trees  retained  is  needed 
as  an  independent  variable,  but  is  itself  an  un- 
known. Successive  percentages  of  trees  are  there- 
fore tested  until  d.b.h.  after  cutting,  number  of 
trees  retained,  and  residual  basal  area  agree  with 
the  growing  stock  level  selected  by  the  program 
user. 

Two  major  loops  in  the  subroutine  provide 
for  the  increasing  and  constant  segments  of  each 
growing  stock  level  curve  (fig.  2).  Limiting  d.b.h. 
for  selection  of  loops  is  10.0  inches  minus  the 
average  change  expected  with  light  partial  cut- 
tings in  stands  just  under  10.0  inches  d.b.h. 

Computed  increases  in  average  d.b.h.  of  dis- 
eased stands  are  half  the  increases  of  similar, 
healthy  stands,  as  occurs  in  actual  dwarf  mistletoe- 
infested  stands,  because  more  of  the  larger  trees 
are  removed  than  would  otherwise  be  the  case. 


Subroutine  SWCUT2 

SWCUT2  computes  the  reduced  average  d.b.h. 
that  results  when  thinning  from  above  in  dwarf 
mistletoe-infested   stands.   The  computations, 
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Figure  2.— Basal  area  after  thinning  in  relation  to 
average  stand  diameter  for  standard  levels  of  grow- 
ing stock. 

which  resemble  SWCUT1,  determine  an  estimated 
post-thinning  d.b.h.  and  compare  the  equivalent 
basal  area  with  the  thinning  level  desired.  Pro- 
gram SWYLD2  duplicates  actual  practice  in  dis- 
eased stands  by  restricting  thinning  from  above 
to  the  initial  entry.  The  program  user  also  has  the 
option  of  requiring  that  all  thinnings  be  from 
below. 


Program  Relationships 

Field  and  office  procedures  used  to  obtain  the 
relationships  described  below  were  similar  to  those 
outlined  previously  for  healthy  (Myers  1971)  and 
dwarf  mistletoe-infested  stands  (Myers  and  others 
1971,  Myers  and  others  1972).  Most  relationships 
in  SWYLD2  are  the  same  as  already  reported 
(Myers  and  others  1972,  Myers  1974).  Additional 
information  on  spread  of  dwarf  mistletoe  from 
overstory  to  understory  was  obtained  from  tem- 
porary plots  in  Arizona  and  New  Mexico. 

Basal  area  and  other  per-acre  values,  average 
stand  diameter,  average  height,  and  site  index  are 
used  as  dependent  and  independent  variables  to 
obtain  the  prediction  equations  used.  The  relation- 
ships shown  in  the  program  listing  (appendix  1) 
contain  only  significant  independent  variables. 

Variable  and  subroutine  names  are  given  in 
capital  letters  in  this  and  the  following  sections. 
Names  of  subscripted  variables  included  in  text 
material  usually  are  not  accompanied  by  symbols 
for  subscripts.  Variables  are  defined  in  the  listing 
of  the  Order  and  Contents  of  the  Data  Deck  and 
in  the  program  listing  (appendix  1). 


Program  statements  of  species-specific  rela- 
tionships must  be  replaced  to  adapt  SWYLD2  to 
other  species  or  conditions.  Suggested  program 
modifications  are  given  in  the  section  headed  Local 
Modification. 


1.  Stocking  After  Cutting 

Stand  density  to  be  left  after  each  partial  cut 
in  healthy  stands,  or  in  diseased  stands  thinned 
from  below,  is  expressed  as  a  relationship  between 
d.b.h.  and  basal  area  (table  1).  Tabular  values, 
computed  as  SQFT  in  subroutine  SWCUT1, 
represent  one  possible  series  of  densities  that  could 
be  used  to  guide  successive  thinnings. 

In  table  1,  basal  area  increases  with  diameter 
until  10.0  inches  diameter  is  reached,  and  remains 
constant  thereafter.  The  designation  "growing 
stock  level  80"  indicates  that  basal  area  is  80.0  ft2 
when  diameter  is  10.0  inches  or  larger,  regardless 
of  what  basal  area  may  be  at  lower  average 
diameters. 

Desired  stand  density  will  vary  with  the 
objectives  of  management,  and  a  family  of  basal 
area-diameter  relationships  is  needed  (see  fig.  2). 
Other  growing  stock  levels  are  named  the  same 
way  as  level  80.  For  example,  level  100  means  that 
reserve  basal  area  will  be  100  ft2  when  d.b.h.  is 
10.0  inches  or  larger.  Basal  areas  for  any  growing 
stock  level  are  computed  by  multiplying  the  values 
for  level  80  in  table  1  by  the  ratio  level/80.  For 
example,  basal  areas  for  level  100  and  diameters 
smaller  than  10.0  inches  are  each  100/80  times  the 
corresponding  basal  areas  of  table  1. 

Desired  levels  are  designated  by  assigning 
values  to  variables  DSTY  and  THIN  on  data  card 
type  4,  as  shown  in  the  listing  of  Order  and  Con- 
tents of  the  Data  Deck.  Each  assigned  value  is  a 
growing  stock  level,  the  basal  area  to  be  left  when 
d.b.h.  after  thinning  is  10.0  inches  or  larger. 

Equations  for  DBHP  in  subroutine  SWCUT1 
also  describe  the  tabulated  values.  In  this  case, 
average  diameters  less  than  10.0  inches  are  esti- 
mated when  basal  areas  and  the  desired  stocking 
level  are  known.  Variables  BREAK  and  BUST 
indicate  points  where  the  relationship  of  diameter 
to  basal  area  has  been  broken  into  segments  for 
convenience  in  regression  analysis. 

Three  equations  for  GSLEVL  in  CUTS, 
SWCUT2,  and  PROJ  are  obtained  from  the 
values  of  table  1,  adjusted  for  various  growing 
stock  levels.  They  compute  the  equivalent  grow- 
ing stock  level  when  basal  area  and  average  diam- 
eter are  known.  The  first  equation  of  each  set  is  a 
linear  extrapolation  to  zero  of  the  data  in  table  1. 


Table  I.  —  Basal  areas  after  partial  cutting  in  relation  to  average  stand  diameter, 

growing  stock  level  80 


Average  Basal 

stand  d.b.h.  area 

after  cutting  per 

(inches)  acre 


Average  Basal 

stand  d.b.h.  area 

after  cutting  per 

(inches)  acre 


Average  Basal 

stand  d.b.h.  area 

after  cutting  per 

(Inches)  acre 


Average  Basal 

stand  d.b.h.  area 

after  cutting  per 

(Inches)  acre 


2.0 
2.1 
2.2 
2.3 
2.4 

2.5 
2.6 
2.7 
2.8 
2.9 


ft2 

12.1 
13.2 
14.4 
15.5 
16.7 

17.9 
19.0 
20.2 
21.3 
22.5 

23.7 
24.8 
26.0 
27.1 
28.3 

29-5 
30.6 
31.8 
32.9 
34.1 


4.0 
4.1 
4.2 
4.3 
4.4 


4.9 

5.0 
5.1 
5.2 
5.3 
5.4 


ft2 

35.2 
36.4 
37.6 
38.7 
39.9 

41 .0 
42.2 
43.4 
44.5 
45.7 

46.8 
47.8 
48.8 
49.8 
50.8 

51.8 
52.8 
53.8 
54.7 
55.7 


6.0 
6.1 
6.2 
6.3 
6.4 

6.5 
6.6 
6.7 
6.8 
6.9 

7.0 
7.1 
7.2 
7.3 
7.4 

7.5 
7.6 
7.7 
7.8 
7-9 


ft2 

56.6 
57.6 
58.5 
59.4 
60.3 

61.2 
62.1 
62.9 
63.8 
64.6 

65.4 
66.2 
67.0 
67.7 
68.5 

69.2 
69.9 
70.6 
71.2 
71.9 


8.0 
8.1 
8.2 
8.3 
8.4 

8.5 
8.6 
8.7 
8.8 
8.9 

9.0 
9.1 
9.2 
9.3 
9.4 

9.5 
9.6 
9.7 
9.8 
9.9 
10.0+ 


ft2 

72.5 
73.1 
73.7 
74.3 
74.8 

75.3 
75.8 
3 
7 
1 


76. 
76. 
77. 


77-5 
77.9 
78.2 
78.5 
78.8 

79.1 
79-3 
79.5 
79.7 
79.8 
80.0 


2.  Diameter  Increase  from  Growth 

Equations  that  predict  future  average  d.b.h. 
of  healthy  stands  (DBHO)  are  obtained  from  plot 
data  that  include  past  and  present  average  diam- 
eters and  other  stand  measures.  Separate  equa- 
tions estimate  growth  at  stand  densities:  (1)  some- 
what less  than  or  within  the  possible  range  of 
managed  stands,  and  (2)  greater  than  any  goals 
of  management. 

Data  from  diseased  plots  provide  the  equa- 
tion for  the  reduced  diameter  growth  caused  by 
dwarf  mistletoe  (DINC).  Plot  values  at  the  begin- 
ning of  the  projection  period  and  the  DBHO 
equation  for  healthy  stands  give  potential  periodic 
growth.  Actual  growth  and  dwarf  mistletoe  rat- 
ings of  the  plots  are  then  used  to  express  actual 
growth  as  a  percentage  of  the  potential  for  various 
degrees  of  infestation. 

Data  for  diameter  growth  equations  may  be 
obtained  from  either  permanent  or  temporary 
plots.  On  permanent  plots,  measurements  must 
be  made  at  intervals  equal  to  the  projection  period 
desired.  Past  diameters  are  obtainable  from  tem- 


porary plot  data  when  present  diameters,  radial 
wood  growth,  and  adjustment  for  bark  growth 
are  known. 

Average  stand  diameter  in  SWYLD2  is  the 
diameter  of  the  tree  of  average  basal  area. 


3.  Diameter  Change  from  Thinning 

Two  sets  of  equations  are  used  to  determine 
average  diameters  after  partial  cutting.  Equations 
for  DBHE  in  subroutine  SWCUTl  estimate  the 
increased  averages  after  cuts  that  emphasize 
removal  of  the  smaller  trees.  The  computed  in- 
crease in  infested  stands  is  half  that  in  healthy 
stands,  to  account  for  the  larger  trees  removed 
in  dwarf  mistletoe  control.  Equations  for  DBHE 
in  subroutine  SWCUT2  compute  the  smaller 
average  diameters  produced  by  thinning  from 
above.  The  latter  situation,  with  emphasis  on 
removal  of  the  larger  trees,  is  the  usual  initial 
thinning  in  infested  stands  if  the  objective  is 
reduction  of  dwarf  mistletoe. 
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Equations  for  post-thinning  diameters  are 
derived  from  data  obtained  in  a  variety  of  stands 
marked  for  thinning.  Thinnings  simulated  on  a 
computer  also  provide  data  for  the  equations 
(Myers  1971). 


4.  Height  Increase  from  Growth 

Equations  in  subroutines  FIRST  and  PROJ 
compute  average  heights  of  dominant  and  codom- 
inant  trees  in  healthy  stands.  Equations  for  HTSO 
and  TEMH  apply  to  stands  of  basal  area  less  than 
or  within  the  range  of  possible  management  goals. 
The  equation  for  HTNEW  applies  to  stands  of 
greater  density.  The  equations  for  HTSO  in 
FIRST  are  bypassed  if  actual  initial  heights  are 
read  in. 

Two  equations  for  reduction  in  height  growth 
due  to  dwarf  mistletoe  are  based  on  data  from 
infested  stands.  The  equation  for  HTPCT  in 
FIRST  computes  any  reduction  during  the  period 
from  stand  origin  to  AGEO  and  is  not  used  if 
initial  heights  are  read  in.  The  equation  for  PCT 
in  PROJ  provides  the  reduction  for  a  single  pro- 
jection period.  Procedures  for  measurement  and 
regression  analysis  of  field  data  parallel  those  for 
the  diameter  growth  equations  described  above. 

Stand  ages  used  in  the  equations  are  averages 
of  dominant  and  codominant  trees.  For  a  two- 
storied  stand,  age  of  the  understory  is  obtained 
from  trees  judged  to  be  potential  dominants  and 
codominants. 


5.  Height  Change  from  Thinning 

Two  equations  for  ADDHT  in  subroutine 
CUTS  compute  the  change  in  average  height  due 
to  partial  cutting.  One  equation  determines  the 
increase  resulting  from  thinning  primarily  from 
below.  The  other  estimates  the  decrease  in  aver- 
age height  that  occurs  when  cutting  from  above 
with  emphasis  on  removal  of  the  larger  trees.  At 
each  cutting,  the  current  value  of  ADDHT  is 
added  algebraically  to  height  before  thinning  to 
obtain  average  height  after  thinning.  Computed 
values  of  ADDHT  are  small  when  thinning  from 
below  because  average  heights  and  their  changes 
refer  to  dominant  and  codominant  trees  only. 


6.  Noncatastrophic  Mortality 

Equations  to  estimate  mortality  in  healthy 
(OUT)  and  diseased  (DIE)  stands  require  values 
for  several  easily  measured  stand  variables.  Data 


for  regression  analyses  came  from  permanent  and 
temporary  plots,  with  and  without  dwarf  mistle- 
toe. As  used  here,  the  term  "noncatastrophic" 
refers  to  gradual  losses  from  competition,  endemic  ' 
insect  populations,  and  dwarf  mistletoe.  Effects 
of  forest  fires,  insect  epidemics,  and  other  abrupt  ; 
changes  are  not  computed  by  these  equations. 

Reduction  in  numbers  of  trees  is  important 
in  unthinned  stands  of  ponderosa  pine,  but  is 
minor  and  erratic  in  healthy,  thinned  stands.  A  3 
prediction  equation  for  mortality  could  not  be 
computed  for  healthy  thinned  stands  with  an  . 
average  d.b.h.  of  10.0  inches  or  larger.  Subroutine 
PROJ  therefore  computes  mortality  reductions  ' 
only  for  mistletoe-free  stands  that  average  less 
than  10.0  inches  d.b.h.  Losses  are  computed  for  f1 
all  infested  stands  with  a  DMR  large  enough  to  ' 
influence  rate  of  mortality.  In  all  computations, 
mortality  is  a  percentage  expressed  as  a  decimal. 

;  lie 

7.  Dwarf  Mistletoe  Infestation 

Ira 

Several  prediction  equations  compute  initial  f 
dwarf  mistletoe  ratings  (DMR)  and  subsequent 
annual  or  periodic  increases.  Data  used  to  obtain  01 
these  equations  come  from  permanent  and  tem-  " 
porary  plots.  Three  ways  of  reporting  initial  degree 
of  infestation  are  provided.  One  alternative  is  to 
enter  DMR  as  part  of  the  input.  Otherwise,  stand    ' ' 
age  at  time  of  initial  infestation  or  the  percentage 
of  trees  infected  may  be  used. 

Initial  ratings  are  computed  by  subroutines  liJ" 
FIRST  or  PROJ,  if  an  initial  DMR  is  not  reported.  * 
The  first  equation  for  an  initial  rating  computes  1  "P 
DMR  from  the  percentage  of  infected  trees  (PINF)  <>P 
in  the  story  or  even-aged  stand.  This  requires  that 
PINF  be  read  in  as  a  second  alternative  way  of  m\ 
reporting  the  degree  of  infestation.  Methods  for  ■ 
determining  PINF  in  southwestern  ponderosa  IP' 
pine  stands  were  described  by  Walters  (1975).  The  1  ^ 
other  three  computations  of  initial  DMR  are  based 
on  another  alternative,  use  of  the  variable  START.  ] 
This  is  the  average  tree  age  when  the  story  or  nine 
even-aged  stand  is  first  infected.  It  is  the  mean  1 
age  when  infection  began  at  various  points  M\ 
throughout  the  story  or  stand.  It  is  not  the  age  lata 
when  earliest  infection  occurred.  Three  equations  ;pprc 
involving  TIME  (derived  from  START)  areneeded,  fecr 
one  for  an  overstory  or  single-storied  stand,  one 
for  understories  with  less  than  3,000  trees  per 
acre,  and  another  for  understories  of  greater  I  St: 
density. 

A  value  of  PINF  is  computed  in  CUTS,  even  \ 
when  other  alternative  measures  of  degree  of  bark 
infestation  have  been  provided.  This  value  for  Ion 
PINF  is  then  used  to  compute  the  buildup  of  tions 
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infestation  over  a  20-year  period  after  DMR  has 
apparently  been  reduced  to  zero  by  partial  cutting. 

Values  of  DMR  computed  from  START  and 
TIME  apply  to  infested  stands  that  have  never 
been  partially  cut.  To  simulate  future  growth  in 
actual  stands  that  have  been  partially  cut,  either 
DMR  or  PINF  should  be  used  to  describe  the 
infestation. 

Four  equations  for  YDM  in  subroutine  PROJ 
provide  the  annual  increase  in  dwarf  mistletoe 
rating  that  is  then  used  to  obtain  the  new  rating 
at  the  end  of  the  projection  period.  The  first  equa- 
tion for  YDM  applies  to  the  increase  in  even-aged 
stands  or  to  overstories.  The  third  and  fourth 
equations  apply  to  understories  and  include  the 
additional  contribution  of  dwarf  mistletoe  seeds 
from  the  overstory.  The  second  equation,  YDM  = 
0.08,  provides  the  constant  increase  in  stands  of 
15  or  less  trees  per  acre.  Spread  of  infection  in 
such  open  stands  is  within-tree  and  not  between 
trees. 

The  two  equations  for  DMR  in  PROJ  that  use 
computed  values  of  YDM  contain  no  species- 
specific  regression  coefficients.  The  weights  of 
I  0.5  and  1.0  assigned  the  periodic  increases  may, 
I  however,  be  different  for  other  combinations  of 
i  host  and  parasite. 


8.  Decrease  in  DMR  from  Thinning 

Two  equations  in  CUTS  compute  the  reduc- 
tion in  dwarf  mistletoe  rating  produced  by  partial 
cutting.  An  initial  cut  in  an  infested  stand  will 
emphasize  removal  of  the  larger  trees  unless  the 
option  is  overridden  by  a  non-zero  value  of  IOPT 
on  data  card  type  4.  The  equation  for  DELDMR 
computes  the  reduction  if  the  initial  thinning  is 
from  above.  The  next  statement  in  the  subroutine 
applies  to  thinnings  from  below,  and  indicates 
that  the  reduction  in  DMR  is  2/3  of  that  possible 
when  thinning  from  above. 

The  equations  were  computed  from  data  ob- 
tained by  simulation  of  thinnings  in  actual  stands. 

Thinnings  are  not  performed  if  the  current 
DMR  exceeds  the  value  of  DMLEV  read  in  by 
data  card  type  4.  A  DMLEV  value  of  3.0  is  usually 
appropriate  for  the  host-parasite  combination 
described  here. 


9.  Stand  Volume  Equation 

Volumes  per  acre  in  total  cubic  feet  inside 
bark  from  groundline  to  tip  (TCF)  are  computed 
from  stand  volume  equations.  Tree  volume  equa- 
tions for  southwestern  ponderosa  pine  (Myers 


1963)  were  used  to  obtain  plot  volumes  for  the 
regression  analysis.  Two  statements  for  TCF 
appear  in  subroutine  SWVOL  because  the  rela- 
tionship is  not  linear  over  the  range  of  D2H  that 
can  appear  in  yield  tables. 

Total  cubic  feet  is  the  only  volume  computed 
directly  by  SWYLD2.  Volumes  in  other  units  are 
obtained  from  total  cubic  volume  and  appropriate 
conversion  factors,  as  explained  below. 


10.  Volume  Conversion  Factors 

Volumes  in  total  cubic  feet  can  be  converted 
to  merchantable  cubic  feet,  board  feet,  square  feet 
of  veneer,  weight  per  acre,  or  other  units.  Appro- 
priate factors  are  used  to  make  conversions.  Sub- 
routine SWVOL  converts  total  cubic  volumes  to 
merchantable  cubic  feet  (with  FCTR)  and  to  board 
feet  Scribner  Rule  (with  PROD).  Footnotes  to  the 
yield  tables  give  minimum  limits  for  tree  d.b.h. 
and  top  diameter. 

Plot  data  used  to  compute  total  cubic  feet 
also  provide  volumes  in  other  units  of  interest 
(Myers  1963).  The  quantity  of  each  unit  per  total 
cubic  foot  is  determined  separately  for  each  plot. 
These  ratios  are  then  used  to  obtain  equations 
for  FCTR  and  PROD.  Three  equations  for  FCTR 
and  two  for  PROD  are  included  in  SWVOL  to 
simplify  regression  analysis. 

Ratios  are  not  computed  for  the  smallest 
average  stand  diameters  because  the  factors  vary 
greatly  in  such  stands.  This  variation  has  no 
important  effect  on  yield  table  construction. 
Merchantable  material  will  not  be  a  part  of  thin- 
nings until  stand  diameter  is  large  enough  for 
conversion  factors  to  be  reliable. 

User-Supplied  Information 

Values  of  38  variables  that  describe  stand 
conditions  or  control  program  execution  are  read 
from  data  cards  supplied  by  the  program  user. 
Twelve  of  these  values  must  be  zero  or  blank,  if 
the  data  apply  to  an  even-aged  stand.  Names, 
input  formats,  sequence  for  punching,  and  defini- 
tions of  these  variables  are  given  in  the  tabulation 
headed  Order  and  Contents  of  the  Data  Deck. 
Twenty  variables  that  require  management  deci- 
sions or  the  collection  of  field  data  are  described 
in  detail  in  this  section. 


1.  Initial  Stand  Measurements 

Values  for  AGEO,  DBHO,  HTSO,  and  DENO 
on  data  card  type  5  are  obtained  in  one  of  two 
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Order  and  Contents  of  the  Data  Deck 


Card     Number  of  Variable 

type         cards  name  Columns  Format 


1 

1 

NSTND 

1-5 

15 

RINT 

6-10 

F5.0 

2 

1 

COMCU 

1-5 

F5.0 

COMBF 

6-10 

F5.0 

3 
4 

1 

per  stand 

1 

per  stand 

NOTE(I) 
ICUT 

1-80 
1-5 

16A5 
15 

DMLEV 

6-10 

F5.1 

DELAY(l) 

11-15 

F5.0 

IOPT(l) 

16-20 

15 

THIN(l) 

21-25 

F5.0 

DSTY(l) 

26-30 

F5.0 

JCYCL(l) 

31-35 

15 

DELAY(2) 

36-40 

F5.0 

IOPT(2) 

41-45 

15 

THIN(2) 

46-50 

F5.0 

DSTY(2) 

51-55 

F5.0 

JCYCL(2) 

56-60 

15 

5 

1 

per  stand 

STND 

SITE 
AGEO(l) 

1-5 

6-10 
11-15 

F5.0 

F5.0 
F5.0 

DBHO(l) 

16-20 

F5.1 

HTSO(l) 

21-25 

F5.1 

DENO(l) 

26-30 

F5.0 

DMR(l) 

31-35 

F5.1 

Description  of  variable 

Number  of  stand  records  to  be  processed  in  a 
single  job. 

Number  of  years  for  which  growth  and  infec- 
tion equations  make  one  projection  of  growth 
or  change.  Value  is  10.0  for  the  equations  in 
appendix  1. 

Minimum  cut  in  merchantable  cubic  feet  to 
be  included  in  total  yields.  Must  be  at  least  1.0. 
Minimum  cut  in  board  feet  to  be  included  in 
total  yields.  Must  be  at  least  1.0. 
Description  of  test  conditions  for  the  stand 
being  processed. 

Control  to  choose:  no  precommercial  thinnings 
(0),  initial  precommercial  thinning  only  (1),  no 
restriction  on  precommercial  thinnings  (2). 
Dwarf  mistletoe  rating  above  which  no  inter- 
mediate cutting  will  be  done. 
Number  of  years  between  initial  age  in  yield 
table  and  tree  age  at  first  cut  in  the  overstory. 
Initial  thinning  in  diseased  overstory  to  be 
from  above  (0)  or  below  (1). 
Growing  stock  level  for  initial  thinning  in  a 
healthy  overstory. 

Growing  stock  level  for  intermediate  cuts 

after  the  first  in  the  overstory. 

Interval  between  intermediate  cuts  in  the 

overstory. 

Number  of  years  between  initial  age  in  yield 
table  and  tree  age  at  first  cut  in  the  under- 
story. 

Initial  thinning  in  diseased  under  story  to  be 
from  above  (0)  or  below  (1). 
Growing  stock  level  for  initial  thinning  in  a 
healthy  understory. 

Growing  stock  level  for  intermediate  cuts 

after  the  first  in  the  understory. 

Interval  between  intermediate  cuts  in  the 

understory. 

Number  of  the  stand. 

Site  index  of  the  stand. 

Age  of  overstory  at  first  entry  in  the  yield 
table. 

Average  d.b.h.  of  the  overstory  at  age 
AGEO(l). 

Average  height  of  overstory  dominants  and 
codominants  at  age  AGEO(l). 
Number  of  live  trees  per  acre  in  the  overstory 
at  age  AGEO(l). 

Average  dwarf  mistletoe  rating  of  overstory 
at  age  AGEO(l). 


1  0 


per  stand 


PINF(l) 

36-40 

F5.1 

START(l) 

41-45 

F5.0 

AGE0(2) 

46-50 

F5.0 

DBH0(2) 

51-55 

F5.1 

HTS0(2) 

56-60 

F5.1 

DEN0(2) 

61-65 

F5.0 

DMR(2) 

66-70 

F5.1 

PINF(2) 

71-75 

F5.1 

START(2) 

76-80 

F5.0 

REGN(l) 

1-5 

F5.0 

VLLV(l) 

6-10 

F5.3 

REGN(2) 

11-15 

F5.0 

VLLV(2) 

16-20 

F5.3 

REGNO) 

21-25 

F5.0 

ways,  depending  on  the  purpose  of  the  computer 
run.  These  stand  descriptors  are  printed  on  the 
first  line  of  the  yield  table;  thus  they  report  initial 
conditions.  Ways  of  obtaining  the  data  are: 

a.  When  the  purpose  is  to  estimate  perform- 
ance of  a  future  or  hypothetical  managed 
stand:  examine  numerous  young  stands 
to  determine  the  average  d.b.h.  and  height 
attained  at  various  ages  for  each  site  class 
and  for  each  of  several  stand  densities.  The 
objective  is  to  describe  stand  conditions 
at  some  selected  age  (AGEO)  if  regenera- 
tion and  subsequent  growth  and  mortality 
progress  according  to  some  goal  set  for  the 
site  class  and  density. 


Percentage  of  overstory  trees  infected  by 
dwarf  mistletoe  at  age  AGEO(l). 
Age  of  overstory  when  dwarf  mistletoe  infec- 
tion begins.  Enter  number  larger  than  largest 
REGN  if  infection  will  not  occur  during  the 
rotation. 

Age  of  understory  at  first  entry  in  the  yield 
table. 

Average  d.b.h.  of  the  understory  at  age 
AGEO(2). 

Average  height  of  potential  dominants  and 
codominants  in  understory  at  age  AGEO(2). 
Number  of  live  trees  per  acre  in  the  under- 
story at  age  AGEO(2). 

Average  dwarf  mistletoe  rating  of  understory 
at  age  AGEO(2). 

Percentage  of  understory  trees  infected  by 
dwarf  mistletoe  at  age  AGEO(2). 
Age  of  understory  when  dwarf  mistletoe  infec- 
tion begins.  Enter  number  larger  than  largest 
REGN  if  infection  will  not  occur  during  the 
rotation. 

Stand  age  when  first  regeneration  cut  will 
occur.  Must  never  be  zero  or  blank.  This  is 
age  for  clearcutting  if  rest  of  card  is  zero  or 
blank. 

Percentage  of  DSTY  to  be  left  atageREGN(l). 
Will  be  zero  with  clearcutting.  Enter  as  a 
decimal. 

Stand  age  at  which  second  regeneration  cut, 
if  any,  will  occur.  Age  for  removal  of  seed 
trees  or  second  cut  of  shelterwood. 
Percentage  of  basal  area  left  after  first  regen- 
eration cut  to  be  left  at  age  REGN(2).  Will  be 
zero  except  for  3-cut  shelterwood.  Enter  as  a 
decimal. 

Stand  age  at  which  third  regeneration  cut,  if 
any,  will  occur.  Final  cut  of  3-cut  shelterwood. 


b.  When  the  purpose  is  to  estimate  future 
performance  of  an  actual  stand:  obtain 
age,  average  d.b.h.,  average  dominant  and 
codominant  height,  and  number  of  trees 
per  acre  by  usual  stand  inventory  pro- 
cedures. 

AGEO  is  the  average  age  of  living  dominant 
and  codominant  trees.  For  an  understory,  age  of 
potential  dominants  and  codominants  is  deter- 
mined separately  from  age  of  the  overstory. 
AGEO  is  the  first  age  that  will  appear  in  the  yield 
table  for  each  stand  and,  unless  a  delay  is  called 
for,  will  be  stand  age  at  time  of  the  first  cutting 
computed  by  SWYLD2.  AGEO  for  actual  stands 
is  age  at  time  of  inventory.  In  the  case  of  hypothet- 
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ical  stands  examined  to  establish  management 
principles,  AGEO  should  usually  be  age  at  time 
of  first  thinning. 

DBHO  is  average  d.b.h.  of  all  live  trees  at  age 
AGEO,  regardless  of  crown  class.  For  two-storied 
stands,  DBHO  is  determined  separately  for  each 
story.  DBHO  should  be  determined  to  the  nearest 
0.1  inch.  Average  d.b.h.  is  the  diameter  of  the  tree 
of  average  basal  area. 

HTSO  is  the  average  height  of  live  dominant 
and  codominant  trees  of  an  even-aged  stand,  or  of 
a  single  story  of  a  two-storied  stand.  As  with 
AGEO,  an  understory  value  is  obtained  from  trees 
judged  to  be  potential  dominants  and  codomi- 
nants.  HTSO  may  be  omitted  from  data  card 
type  5  if  a  hypothetical  stand  with  a  density  below 
growing  stock  level  160  is  being  examined.  Values 
for  HTSO  are  computed  if  zero  or  blank  fields 
appear  on  the  data  cards.  For  actual  stands, 
measured  heights  should  be  used  to  account  for 
the  past  effects  on  height  growth  of  stand  den- 
sity and  disease. 

DENO  is  the  number  of  live  trees  per  acre  in 
an  even-aged  stand  or  in  each  story  of  two-storied 
stands.  All  trees  over  4.5  feet  tall  are  counted, 
regardless  of  crown  class.  These  are  the  same 
trees  used  to  obtain  DBHO. 


2.  Initial  Dwarf  Mistletoe  Conditions 

Three  situations  involving  time  of  initial 
infestation  can  be  examined  for  actual  or  hypo- 
thetical stories  or  stands.  For  each,  a  non-zero 
value  of  one  variable  is  punched  in  data  card  5  and 
the  values  of  the  other  two  disease  variables  are 
zero  or  blank.  The  situations  are: 

1.  The  story  or  stand  is  already  infested  at 
age  AGEO,  the  first  age  in  the  yield  table. 
Either  the  DMR  or  the  percentage  of  trees 
infected  is  known. 

2.  Infestation  begins  at  any  time  before  or 
during  the  period  simulated,  and  story  or 
stand  age  at  this  time  is  known. 

3.  The  story  or  stand  remains  healthy 
throughout  the  period  simulated. 

An  infestation  already  present  at  age  AGEO 
is  described  in  either  of  two  ways: 

1.  DMR— average  dwarf  mistletoe  rating 
(see  fig.  1)  at  age  AGEO  of  all  trees  of  an 
actual  story  or  stand  or  a  selected  value 
assigned  to  a  hypothetical  stand.  A  non- 
zero value  of  DMR  on  the  data  card  causes 
statements  that  compute  initial  DMR  to 
be  bypassed. 


2.  PINF— percentage  of  the  live  trees  infected 
with  dwarf  mistletoe  in  an  actual  or  hypo- 
thetical story  or  stand  at  age  AGEO.  DMR 
is  later  computed  from  PINF. 

An  infestation  that  begins  at  a  known  stand 
or  story  age,  before  or  after  AGEO,  is  described 
by  the  variable  START.  This  is  the  average  age 
of  the  dominant  and  codominant  trees  when  the 
story  or  stand  is  first  infected.  Refer  to  item  7  of 
"Program  Relationships"  for  a  detailed  definition 
of  START.  Preferably,  START  should  be  used 
only  when  simulating  growth  in  hypothetical 
stands.  DMR  is  computed  from  START  only  if 
both  the  DMR  and  PINF  fields  of  data  card  type  5 
are  zero  or  blank. 

All  computations  of  DMR  must  be  bypassed 
if  the  story  or  stand  remains  mistletoe-free 
throughout  the  period  simulated.  This  is  done  by 
entering  zero  values  for  both  DMR  and  PINF  and 
a  value  for  START  greater  than  the  largest  value 
of  the  three  REGN  variables.  For  example,  if  with 
2-cut  shelterwood  REGN(l)  equals  110,  REGN(2) 
equals  130,  and  REGN(3)  equals  0,  then  START 
must  have  some  value  greater  than  130. 


3.  Site  Index 

The  site  index  value  read  in  for  each  stand  is 
based  on  mean  height  of  dominant  and  codominant 
trees  at  a  base  age  of  100  years.  Tables  and  graphs 
published  by  Meyer  (1938)  provide  the  relation- 
ships used  to  determine  each  index. 

Site  index  values  used  in  growth  computa- 
tions for  actual  stands  must  be  reliable  indicators 
of  relative  productivity.  Trees  in  dense  stands  or 
with  any  disease  or  deformity  severe  enough  to 
reduce  height  growth  do  not  furnish  reliable  infor- 
mation. In  such  cases,  site  index  can  be  obtained 
from  suitable  nearby  stands  or  from  soil  variables. 


4.  Intermediate  Cuts 

Seven  variables  on  data  card  type  4  control 
the  execution  of  intermediate  cuts.  Eleven  fields 
are  provided  for  them,  since  four  variables  require 
separate  entries  for  each  story  of  two-storied 
stands.  The  purpose  of  each  variable  is  as  follows: 

1.  ICUT— control  of  the  number  of  noncom- 
mercial thinnings  that  will  be  simulated.  An  input 
value  of  zero  causes  all  noncommercial  thinnings 
to  be  bypassed.  Assign  ICUT  a  value  of  one  if 
only  the  initial  thinning  can  be  noncommercial. 
An  input  value  of  two  results  in  all  scheduled 
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thinnings  being  performed,  regardless  of  how 
many  may  be  noncommercial.  Minimum  volumes 
for  commercial  cuts  are  entered  as  COMBF  and 
COMCU  on  data  card  type  2.  Variables  COMBF 
and  COMCU  must  be  assigned  values  of  at  least 
1.0,  regardless  of  the  value  of  ICUT. 


2.  THIN— growing  stock  level  for  initial 
thinning  in  healthy  stands.  Levels  are  designated 
in  the  manner  shown  in  figure  2.  They  will  fre- 
quently be  80,  90,  or  100  for  southwestern  pon- 
derosa  pine,  but  may  range  from  40  to  160. 


7.  IOPT— control  on  whether  or  not  the  initial 
thinning  in  an  infested  stand  will  be  from  above. 
Thinning  from  above  emphasizes,  but  is  not 
restricted  to,  removal  of  the  larger  trees.  Such 
cutting  is  usual  practice  in  infested  stands  because 
the  larger  trees  usually  support  the  greater 
amounts  of  dwarf  mistletoe.  IOPT  is  assigned  a 
value  of  zero  if  initial  thinning  is  to  be  from  above 
and  a  value  greater  than  zero  if  the  initial  thinning 
is  to  be  from  below.  Regardless  of  the  value  of 
IOPT,  thinning  will  not  be  performed  unless  all 
restrictions  on  thinning  are  met. 


3.  DMLEV— maximum  dwarf  mistletoe  rat- 
ing above  which  intermediate  cuttings  will  not  be 
performed.  A  value  of  3.0  will  often  be  used  for 
ponderosa  pine.  Experience  indicates  that  thin- 
nings cannot  materially  reduce  effects  of  the 
infestation  or  improve  stand  growth  when  DMR 
exceeds  3.0. 


4.  JCYCL— interval  between  intermediate 
cuts,  as  punched  on  card  type  4.  JCYCL  is  later 
redefined  in  subroutine  CUTS  when  an  interval 
between  regeneration  cuts  is  needed.  Values 
assigned  to  JCYCL  must  be  a  multiple  of  the 
projection  period  of  the  growth  and  mortality 
equations  (RINT).  For  example,  the  equations  of 
SWYLD2  project  d.b.h.  and  other  measures  for 
10-year  periods,  and  RINT  has  a  value  of  10.0. 
The  value  of  JCYCL  on  card  type  4  can  therefore 
be  10,  20,  30,  etc.  years. 


5.  DSTY— growing  stock  level  for  interme- 
diate cuts  after  the  first  cut.  Discussion  of  choice 
of  level  for  THIN,  above,  also  applies  to  DSTY. 
Levels  assigned  to  THIN  and  DSTY  may  be  equal. 

6.  DELAY— number  of  years  between  the 
first  age  reported  in  the  yield  table  (AGEO)  and 
age  of  the  story  or  stand  at  time  of  the  first  simu- 
lated cutting.  The  value  of  DELAY  must  be  a 
multiple  of  the  value  of  RINT.  The  variable 
DELAY  schedules  the  first  cut  but  does  not 
assure  that  it  actually  will  be  performed.  Restric- 
tion on  noncommercial  thinning,  basal  area  less 
than  the  growing  stock  goal,  or  a  value  of  DMR 
larger  than  DMLEV  will  prevent  or  postpone 
thinning  regardless  of  the  value  of  DELAY. 
Initial  cutting  will  take  place  at  age  AGEO  if 
DELAY  is  assigned  a  value  of  zero  and  if  other 
restrictions  on  cutting  do  not  apply. 


5.  Regeneration  Cuts 

Entries  on  data  card  type  6  determine  the 
method  of  regeneration  to  be  simulated,  stand 
age  when  each  cut  will  be  performed,  and  the 
basal  area  to  be  left  for  seed  trees  or  shelterwood. 
Stand  ages  at  regeneration  cuts  must  equal  the 
sum  of  AGEO  plus  appropriate  multiples  of  the 
value  of  RINT.  Regardless  of  regeneration  method, 
REGN(l)  will  never  be  assigned  a  value  of  zero 
on  data  card  type  6. 

For  clearcutting,  stand  age  at  time  of  final 
cutting,  or  REGN(l),  is  the  only  entry  needed  on 
the  card. 

Seed  tree  cutting  requires  that  values  for 
REGN(l),  VLLV(l),  and  REGN(2)  be  punched  in 
card  type  6.  REGN(l)  is  the  stand  age  at  first 
regeneration  cutting,  and  REGN(2)  is  stand  age 
when  the  seed  trees  are  removed.  The  interval 
between  AGEO  and  REGN(l)  and  the  interval 
between  REGN(l)  and  REGN(2)  must  each  be  a 
multiple  of  RINT.  VLLV(l)  is  the  percentage  of 
the  growing  stock  goal,  DSTY,  to  be  left  for  seed 
trees.  VLLV(l)  is  entered  as  a  decimal. 

Shelterwood  cuttings  are  controlled  in  the 
same  manner  as  seed  tree  cuts  except  that  up  to 
three  regeneration  cuts  are  possible.  The  restric- 
tions on  formats  and  intervals  described  above 
also  apply  here.  Two-cut  shelterwood  requires 
that  values  be  assigned  to  the  same  three  variables 
used  with  the  seed  tree  method,  described  above. 
For  three-cut  shelterwood,  REGN(l)  and  REGN(2) 
are  stand  ages  at  the  two  removal  cuts  and 
REGN(3)  is  stand  age  at  the  final  cut.  VLLV(2) 
is  the  percentage  of  the  shelterwood  basal  area 
left  at  age  REGN(l)  that  will  be  retained  at  age 
REGN(2).  For  example,  assume  that  DSTY  is 
100  square  feet  and  that  VLLV(l)  and  VLLV(2) 
are  0.5  and  0.4,  respectively.  Basal  area  of  the 
shelterwood  at  age  REGN(l)  will  be  50  square 
feet.  Residual  basal  area  at  age  REGN(2)  will  be 
20  square  feet. 
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Local  Modifications 

This  section  may  be  skipped  by  readers  who 
are  not  modifying  SWYLD2  for  local  use.  If  the 
adaptation  is  to  apply  to  healthy,  single-storied 
stands  only,  it  will  be  more  efficient  to  modify 
program  PONYLD  (Myers  1971)  rather  than 
SWYLD2. 

Complete  modification  of  SWYLD2  for  other 
host-dwarf  mistletoe  combinations,  regions,  thin- 
ning specifications,  and  utilization  standards 
requires  replacement  of  65  statements  that  com- 
pute values  of  24  variables.  The  statements  are 
described  in  the  section  headed  "Program  Rela- 
tionships." The  65  replacements  require  that 
coefficients  for  46  different  relationships  be  com- 
puted. The  difference  of  19  is  due  to  duplication 
among  the  65  replaceable  statements.  The  46  dif- 
ferent relationships  include  several  cases  where, 
for  simplicity  of  regression  analysis,  a  single 
variable  is  computed  by  two  or  three  separate 
relationships. 

The  65  replaceable  statements  described  below 
do  not  include  the  IF  statements  needed  to  select 
one  appropriate  statement  from  among  the  two 
or  three  provided  for  some  variables.  The  number 
of  IF  statements  to  be  modified  or  removed  will 
depend  on  the  results  of  each  regression  analysis. 

Statements  to  be  replaced  for  complete  modi- 
fication are  in  six  subroutines,  as  follows: 

FIRST  has  seven  replaceable  statements: 
four  for  DMR,  two  for  HTSO,  and  one  for  HTPCT. 

CUTS  has  eight  replaceable  statements:  three 
for  GSLEVL,  two  for  DELDMR,  two  for  ADDHT, 
and  one  for  PINF. 

PROJ  has  23  replaceable  statements: 

Variable     Statements     Variable  Statements 


DMR 

6 

DIE 

1 

DENT 

1 

OUT 

1 

YDM 

4 

TEMH 

2 

GSLEVL 

3 

HTNEW 

1 

DBHO 

2 

PCT 

1 

DINC 

1 

The  first  four  statements  for  DMR  contain 
the  same  coefficients  as  the  four  in  subroutine 
FIRST.  The  last  two  statements  for  DMR  have 
weighting  factors  of  0.5  and  1.0  that  could  vary 
by  species.  The  IF  statement  using  DENT, 
immediately  following  the  statement  labeled  20, 
contains  a  constant  that  may  be  species-specific. 
Three  statements  for  YDM  resemble  three  of  the 
statements  for  DMR,  but  do  not  include  the 
variable  TIME.  The  other  statement  for  YDM 
sets  a  constant  value  of  0.08  that  may  be  species- 
specific.  Statements  for  GSLEVL  are  the  same 


in  PROJ  and  CUTS.  Two  statements  for  TEMH 
have  the  same  coefficients  as  the  statements  for 
HTSO  in  subroutine  FIRST.  The  remaining  seven 
replaceable  statements  in  PROJ  are  not  dupli-  j 
cated  in  PROJ  or  in  other  subroutines. 

SWVOL  has  seven  replaceable  statements:  t 
two  for  TCF,  three  for  FCTR,  and  two  for  PROD. 
None  of  them  are  duplicated  elsewhere  in  the 
program. 

SWCUT1  has  14  replaceable  statements,  none  % 
of  which  are  duplicated  elsewhere.  Long  state- 
ments for  DBHE  and  PDBHE  (log-log  relation- 
ship for  DBHE)  each  appear  twice.  Two  short!  | 
identical  statements  for  DBHE  just  before  the  i 
statements  labeled  20  and  70,  contain  a  constant  :Ii 
0.5  that  may  be  species-specific.  With  other 
host-parasite  combinations,  thinning  from  below 
in  diseased  stands  may  not  produce  half  the  | 
change  in  d.b.h.  that  thinning  from  below  produces  „, 
in  healthy  stands.  The  remaining  eight  replace- J  j,f 
able  statements  are  not  duplicated  in  SWCUTl  ,3 
or  elsewhere.  They  include  one  statement  each 
for  BREAK  and  BUST,  three  statements  fort  j,, 
SQFT,  and  three  statements  for  DBHP. 

SWCUT2  has  six  replaceable  statements,  four|  j, 
of  which  are  duplicated  elsewhere.  The  three  equa-i  ^ 
tions  for  GSLEVL  have  the  same  coefficients  as,  | 
those  in  CUTS  and  PROJ.  The  equation  for 
DBHE  just  before  the  statement  labeled  30  also,  ^ 
appears  in  SWCUTl.  The  two  long  statements  .!. 
for  DBHE  (one  as  PDBHE  for  a  log-log  relation- 
ship) compute  the  reduced  d.b.h.  after  thinning 
from  above,  and  do  not  appear  elsewhere. 

Seventeen  statements  relating  to  stocking 
need  not  be  replaced  if  the  method  of  expressing 
stocking  levels  shown  in  table  1  and  figure  2  is 
retained.  These  17  include  nine  statements  for 
GSLEVL  in  CUTS,  SWCUT2,  and  PROJ.  Other 
retainable  statements  are  in  SWCUTl:  three  for 
DBHP,  one  each  for  BREAK  and  BUST,  and  ; 
three  for  SQFT. 

Several  statements  in  subroutine  TABLE 
that  write  table  headings  and  footnotes  must  be 
modified  to  conform  to  changes  in  species,  units  f 
of  measure,  or  utilization  standards. 

i ; 

A  Sample  Problem 

The  problem  described  below  demonstrates 
most  of  the  computations  made  by  SWYLD2  and  j  J 
the  printed  results  obtained.  It  illustrates  some  [ 
of  the  questions  that  may  be  asked  and  the  infor- 
mation that  will  be  provided.  For  brevity,  only 
three  of  the  many  possible  stand  conditions  are  ! 
examined.  The  example  also  serves  as  a  test  1  j 
problem  for  use  in  adapting  the  source  program 
to  locally  available  computing  facilities.  t 
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Assume  a  forest  composed  of  even-aged  and 
two-storied  stands  of  ponderosa  pine  with  differ- 
ing degrees  of  infestation  by  dwarf  mistletoe. 
Some  of  the  many  questions  the  manager  of  such 
a  forest  would  have  and  their  relation  to  manage- 
ment include: 

1.  What  growth  can  be  expected  in  a  hypo- 
thetical even-aged  stand,  without  dwarf  mistle- 
toe, that  is  to  be  managed  according  to  a  selected 
set  of  management  controls?  Initially,  many 
yield  tables  would  be  produced  to  examine  the 
effect  of  variations  in  length  of  cutting  cycle  and 
other  controls  on  yields  and  on  the  number  of 
precommercial  thinnings  required. 

2.  How  would  this  performance  be  affected 
by  a  dwarf  mistletoe  infestation  that  started  when 
the  stand  was  10  years  old?  Tables  of  this  type 
provide  information  on  impacts  of  the  disease, 
free  of  the  accidents  of  past  history  of  actual 
stands. 

3.  What  yields  can  be  expected  from  an  exist- 
ing two-storied  stand  infested  with  dwarf  mistle- 
toe if  the  overstory  is  removed  now?  This  and 
the  next  question  relate  to  determination  of  the 
best  way  of  converting  an  undesirable  stand  con- 
dition to  something  better. 

4.  What  would  be  the  effect  on  yield  if  ques- 
tion 3  were  modified  to  delay  overstory  removal 
for  20  years? 


For  purposes  of  this  example,  assume  that 
the  values  of  several  management  controls  have 
been  standardized  for  the  working  circle.  Cutting 
cycle  length  is  20  years.  Initial  thinnings  in  young 
stands  reduce  density  to  growing  stock  level  120. 
Other  intermediate  cuts  are  to  growing  stock 
level  90.  Only  the  first  thinning  may  be  noncom- 
mercial. Minimum  commercial  volumes  per  acre 
are  320  cubic  feet  to  a  4-inch  top  and  1,500  board 
feet.  Regeneration  is  by  two-cut  shelterwood  with 
a  removal  cut  at  stand  age  110  and  the  final  cut 
at  age  130.  The  manager  expects  that  his  pro- 
cedure for  regeneration  cuts  will  result  in  a  new 
stand  with  950  trees  per  acre  at  age  30.  This  new 
crop  will  average  4.8  inches  in  diameter  on  areas 
of  site  index  70. 

The  data  deck  consists  of  the  18  cards  shown 
in  figure  3.  The  first  two  cards  shown  are  the 
type  1  and  type  2  cards,  respectively,  read  by  the 
main  program.  The  next  four  cards  are  the  type  3, 
4,  5,  and  6  cards  used  to  answer  the  first  of  the 
questions  listed  above.  The  remaining  12  cards 
are  the  type  3,  4,  5,  and  6  cards  for  the  other  three 
questions. 

Yield  tables  produced  by  SWYLD2  to  answer 
the  four  questions  given,  and  with  the  data  deck 
described  above,  are  reproduced  in  appendix  2. 

Yield  tables  can  be  used  in  many  ways  to 
assist  in  decisionmaking.  For  many  purposes, 
yields  of  healthy  stands  will  be  desired  so  that 


COLUMN  NUMBERS 

12345678 
12345678 90 123456789Q123456789012345678901234567890 1234567890 12 34567890 123<+567890 

AR0  TYPE 


1  4  10 

2  320  1500 

3  HYPOTHETICAL  HEALTHY  EVEN-AGED  STANO.  BEGIN  TABLE  AT  AGE  30. 
k  1       30  0         0  120       90       20  0         0          0         0  0 

5  1        70  30       46  0     950         0  0     999         0  0         0          0         0         0  0 

6  110     500  130         0  0 

3  HYPOTHETICAL  EVEN-AGED  STANO.    INFESTATION  STARTS  AT   AGE  10. 

4  1       30  0         0  120       90       20  0         0          0         0  0 

5  2       70  30       48  0     950         0  0       10          000000  0 

6  110     500  130  0  0 

3  TWO-STORIED  INFESTED  STAND    HITH   IMMEDIATE  OVERSTORY  REMOVAL. 

4  1        30  0  0  120       90       20  0         0      120       90  20 

5  3       70  150     161  750       20       48  0         0       30       41     210     600       10         0  0 

6  110     500  130  0  0 

3  TWO-STORIED  INFESTED  STANO.   OVERSTORY  REMOVAL   DELAYEO  20  YEARS. 

4  1       30  20         0  120       90       20  20         0     120       90  20 

5  i»       70  150     161  750       20       48  0         0        30       41     210     600       10         0  0 

6  110     500  130         0  0 


Figure  3.— Data  deck  for  sample  problem. 
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long-range  goals  can  be  determined.  Yields,  num- 
ber of  noncommercial  cuts,  number  of  scheduled 
cuts  that  cannot  be  made,  and  size  of  the  average 
tree  are  some  of  the  values  produced.  Money 
yields  and  rates  earned  can  be  computed  if  neces- 
sary data  on  costs  and  stumpage  values  are  avail- 
able. Stand  ages  at  culmination  of  mean  annual 
increment,  and  rates  earned  can  help  the  manager 
determine  suitable  rotations  for  his  working 
groups. 

With  a  disease  as  important  as  dwarf  mistle- 
toe, it  is  also  necessary  to  project  the  growth  of 
diseased  stands.  Two  important  types  of  informa- 
tion can  be  obtained  from  yield  tables  that  include 
disease  effects.  A  comparison  of  projected  yields 
of  healthy  and  diseased  stands  indicates  the 
impact  of  the  disease  on  the  timber  resource. 
Yield  tables  for  actual  diseased  stands  indicate 
the  types  of  treatment  that  may  be  practical,  and 
give  an  estimate  of  the  expected  benefits. 
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APPENDIX  1 

Listing  of  Program  SWYLD2 


TO  COMPUTE  ANO  PRINT  YIELD  TABLES  FOR  EVEN-AGED  ANO  TWO-STORIED  STANDS 
OF   SOUTHWESTERN  PONQEROSA   PINE  WITH    OR   WITHOUT    DWARF  HISTLETOE. 

DEFINITIONS  OF  VARIABLES. 
SUBSCRIPTS  - 

1=1  FOR  OVERSTORY,  1=2  FOR  UNOERSTORT. 
J  =  lt  Z%  OR  3  TO  NUMBER  REGENERATION  CUT. 
K   =    AGE  CLASS  AT    T I  HE    OF  OPERATION. 

ADDHT (I)    =   INCREASE   OR   DECREASE   IN  AVERAGE  STAND   HEIGHT  FROM 

PARTIAL  CUTTING. 
AGEOU)    =    INITIAL   AND    SUBSEQUENT   AGES    IN   YIELD  TABLE. 
BASO(I)    =    BASAL    AREA    PER   ACRE    BEFORE    PARTIAL  CUTTING. 
BAST(I)    =    BASAL  AREA  PER  ACRE   AFTER   PARITAL  CUTTING. 
BREAK    =  POINT   WHERE   GROWING   STOCK  RELATIONSHIP   DF    01 A METER  ON 

BASAL  AREA   IS  BROKEN  INTO  SEGMENTS. 
BUST   =   POINT    WHERE    GROWING  STOCK  RELATIONSHIP   OF   OIAHETER  ON 

BASAL   AREA   IS   BROKEN  INTO  SEGMENTS. 
CBAS ( 1 1 K)    =  BASAL   AREA  REMOVE  D • 
CBF ( I  * K )    =  H  80.  FT.  REMOVED. 
COENd.K)    =  NUMBER   OF    TREES  REMOVEO. 

CHNG  =   PERIOOIC  INCREASE   IN  AVERAGE   DOH  AND  CODDN  HEIGHT. 

CMCF(I»<)    =  KERCH.   CU.   FT.  REMOVED. 

COMBF    =   MINIMUM  COMMERCIAL  CJT,    BOA  RO  FEET. 

COMCU   =    MINIMUM  COMMERCIAL  CUT,    CU.  FT. 

CTOT  ( I « K  >    =  TOTAL   CU.   FT.  REMOVEO. 

DBHE    =    TEMPORARY   VARIABLE   FOR    AVERAGE    ST ANO  O.3.H.    AFTER  PART  I AL 
CUTTING. 

DBHO(I)    =    AVERAGE   STANO  O.B.H.    BEFORE   PARTIAL   C J T TING . 

OBHP   =   AVERAGE  O.B.H.    FROM  GROWING  STOCK  LEVEL  CURVES. 

DBHT(I)    =    AVERAGE   STANO  O.B.H.    AFTER   PARTIAL  CUTTING. 

OELA  Y ( I )    =   NUMBER  OF    YEARS  BETWEEN    INITIAL  AGE    AND   FIRST  PARTIAL 

CUTTING,    IF   FIRST    PARTIAL   CUTTING    IS  NOT    OONE    AT  ONCE. 
OELOMR    =    CHANGE    IN   D.M.R.    DUE    TO   PARTIAL  CUTTING. 
OENO(I)    =    TREES  PER  ACRE   BEFORE    PARTIAL  CUTTING. 
DENT  1 1 )    =    TREES  PER  ACRE   AFTER  PARTIAL  CUTTING. 

DIE  =    TREES   LOST    IN   DISEASED   STANDS    IN  RINT    YEARS,    IN  PERCENT. 
OINC   =   PER  IODIC   INCREASE   IS  AVERAGE   O.B.H.    IN  STANDS  WITH  D.M.R. 

GREATER  THAN  3.5. 
DLEVlIt    =    GROWING   STOCK  LEVEL    FOR   INTERMEDIATE   CUTS   AFTER  THE 

FIRST  AND  FOR  REGENERATION  CUTS  BY  SEEO  TREE  OR  SHEL  TERWOOQ . 
OMLEV  =  O.H.R.  ABOVE  WHICH  NO  INTERMEDIATE  CUTTING  WILL  BE  OONE. 
OMR ( I )    =  DWARF   MISTLETOE    INFECTION  RATING. 

DMR  T ( I  I    =    MAXIMUM    INFECTION  EXPECTED    IN   STANOS   AFTER  PARTIAL 
CUTTING. 

OSTY(I)    =   GROWING  STOCK  LEVEL   FOR  INTERMEDIATE   3UTS   AFTER  THE 
FIRST. 

FCTR  =   FACTOR  TO  CONVERT  TOTAL   CU.   FT.    TO  MERCH.    CU.  FT. 
GSLEVL   =   CURRENT   ACTUAL  GROWING  STOCK  LEVEL. 

HTNEH  =   ESTIMATED  AVERAGE  TREE   HEIGHT   IN  STANOS  WITH  GROWING  STOCK 

LEVEL   GREATER   THAN  160. 
HTPCT    =   PERCENTAGE    OF    POTENTIAL    HEIGHT    GROWTH    ACTUALLY  ATTAINED 

BY    INFESTED  STANDS. 
HTSO(I)    =    AVERAGE    TREE    HEIGHT    BEFORE    PARTIAL  CUTTING. 
HTST (  1 1    =    AVERAGE    TREE    HEIGHT    AFTER    PARTIAL  CUTTING. 
ICUT    =   FLAG   WITH    INPUT    VALUE   DEFINED   BELOW   TO    CONTROL  THINNINGS. 

0  -    NO  NONCOMMERCIAL    THINNINGS  ALLOWED. 

1  -   PRECOMHE  RC I A  L    INITIAL    THINNING   ALLOWED.  NONCOMMERCIAL 

SUBSEQUENT    THINNINGS    NOT    ALL  OWEO. 

Z    -    ALL    THINNINGS  MAY   BE  NONCOMMERCIAL. 
IK(I>    =  NUMBER  OF  GROWTH   PROJECTION  PERIODS  IN  JCYCL(I). 
IOPT(I)    =   FLAG  TO   ALLOW   (0)   OR  PREVENT    (It    INITIAL   THINNING  FROM 

ABOVE    IN  STANDS  WITH   DWARF  MISTLETOE. 
IREGN(I)    =   FLAG    WITH    A   VALUE   OF    ONE    IF   A   REGENERATION  CUT    li  BEING 

PERFORMED. 

ISANSdl    =  FLAG   WITH  A  VALUE  EQUAL   TO  THE  NUMBER   OF  GROWTH 

PROJECTIONS  SINCE  A   CJT   WHICH  REDUCES  O.M.R.    TO   ZERO  IN 

INFESTED  STANOS. 
JCYCHII    =   INTERVAL   BETWEEN  INTERMEDIATE   OR  REGENERATION  CUTS. 
KDEL  <  1 1    =    NUMBER  OF    PROJECTION   PERIODS   BEFORE    INITIAL  PARTIAL 

CUTTING   IS  PERFORMED. 
KOL    =   COUNTER  TO  CONTROL    STORAGE    IN    ARRAYS    PRINTED   IN   YIELO  TABLES 
KOUNT ( 1 1    =    NUMBER  OF  GROWTH   PROJECTIONS   SINCE   LAST  CUT. 
KPER   =    NUMBER  OF    PROJECTION  PERIODS    FROH   INITIAL    AGE    TO  CURRENT 

AGE. 

KSTEP(I)    =   INDICATOR  WITH  VALUE   OF   ONE   IF  CURRENT   PARTIAL  CUTTING 
IS  FROM  BELOW   AND    TWO  IF    CURRENT    PARTIAL    CJTTING    li  FROM 
ABOVE. 

KTRII)    =    INOICATOR   WITH  VALJE   GREATER   THAN  ZERO   IF  A  SCHEOULEO 
THINNING  HAS    BEEN  SKIPPED   BECAUSE   D.M.R.    13    GREATER  THAN 
OHLEV,  STANO   IS   ALREADY  AT   OR  BELOW  SPECIFIEO   STOCKING,  OR 
THINNING  010   NOT    HEET    IMPOSED    MERCHANTABILITY  STANOAROS. 

NOTE (L )    =    DESCRIPTION   OF    TEST   CONOITIONS    FOR  STANO  OUTPUT. 

NSTND   =   NUMBER   OF    STANOS    BEING    PRQCESSEO   BY  RUN. 

OUT    =   PERCENTAGE    MORTALITY    IN    HEALTHY  STANDS. 

PCT    =   PERCENTAGE   OF    POTENTIAL    PERIOOIC    HEIGHT    GROWTH  ACTUALLY 

ATTAINED  BY   INFESTED  STANOS. 
PINFUf    =    PERCENTAGE  OF  TREES   INFECTEO  WITH  OHARF  HISTLETOE. 
PRET   =  PERCENTAGE   OF   TREES  RETAINEQ   AFTER  PARTIAL  CUTTING. 
PROO  =   FACTOR  TO  CONVERT  TOTAL   CU.   FT.    TO  BO.  FT. 
REGN ( J 1    =   STAND  AGE   WHEN  REGENERATION  CUT  J  OCCURS. 
REST (I)    =   GROWING  STOCK  LEVEL   FOR  CURRENT   INTERMEDIATE  CUT. 
RINT    =    NUMBER  OF   YEARS   FOR  WHICH    A    SINGLE   PRO  JE  3  T I  ON   IS  MADE. 
ROTA    =  FINAL   AGE    IN    YIELD  TABLE. 

ROTAOS  =   AGE   OF  OVERSTORY  AT   TIME  OF   FINAL  REMOVAL   IF  OVERSTORY  IS 
TO  BE   RETAINED   PAST    LARGEST  REGN ( J )  . 


S8AS    =    TOTAL    BASAL    AREA   OF   STAND,    INCLUDING  BOTH  STORIES. 
SB0(I1    =    SUM    OF   BOARD   FEET   FROM   ALL    CUTS   WITH    YIELO   OF   COHBF  OR 

LARGER. 
SITE   =   SITE  INOEX. 

SMC(I)    =    SUM    OF    MERCH.    CU.    FT.    FROM    ALL    CUTS   WITH    YIELD    OF  COMCU 
OR  LARGER. 

SQFT    =    BASAL    AREA   AS   A  FUNCTION   OF    AVERAGE    O.B.H.    FOR  GROWING 

STOCK  LEVEL  80. 
ST  ART  1  I )    =   AGE    OF    STORY    AT   TIME    OF    INITIAL  lNF£:TI0N. 
STF(I)    =   SUM  OF   TOTAL   CU.   FT.   FROM  ALL  CUTS. 
STND   =    STAND    IDENTIFICATION  N JMBER. 
T AGE ( I » Kl    =  STANO  AGE   AFTER  CUTTING. 
TBAS ( I  * Kl    =  BASAL   AREA  AFTER  CUTTING. 
TBF(I.K)    =  M  8D.  FT.   AFTER  CUTTING. 
TCF   =    TEMPORARY    VARIABLE    FOR    TOTAL    CU.  FT. 
TOEN(I»K>    =   NUMBER    OF    TREES  AFTER  CUTTING. 
TDM(I.K)    =    AVERAGE    O.B.H.    AFTER  CUTTING. 

TOMR ( I  ,  K)    =   DWARF    HISTLETOE    INFECTION    RATING  AFTER  CUTTING. 
TEHH(L)    =    TEMPROR A Y    VARIABLES    USED   IN   COMPUTATION  OF  PERIOOIC 
HEIGHT  GROWTH. 

THIN(I)    =    GROWING  STOCK  LEVEL    FOR  INITIAL    THINNING,  HEALTHY 
STANDS. 

THT ( I , K )    =    AVERAGE    HEIGHT  AFTER  CUTTING. 

TIME   =   NUMBER  OF   YEARS  STANO   INFESTED   BY   DWARF  MISTLETOE. 

THCF(I,K)    =   MERCH.    CU.   FT.    AFTER  CUTTING. 

TTOT(UK)    =  TOTAL   CU.   FT.   AFTER  CUTTING. 

UAGE(I.K)    =  STANO   AGE   BEFORE  CUTTING. 

UB AS  <  I »  K)    =   BASAL    AREA  BEFORE  CUTTING. 

U3F(I,K>    =   M   BO.   FT.    BEFORE  CUTTING. 

UDEN(I,K>    =  NUMBER   OF    TREES   BEFORE  CUTTING. 

UDM(I.K)    =   AVERAGE   D.B.H.   BEFORE  CJTTING. 

UOMR(I.K)    =   DWARF    HISTLETOE    INFECTION    RATING   8E=0RE  CUTTING. 
UHTtI,K>    =   AVERAGE    HEIGHT   BEFORE  CUTTING. 
UHCF(I,K>    =  MERCH.   CU.   FT.   BEFORE  CUTTING. 
ur0T!I,K>    =  TOTAL  CU.    FT.  BEFORE  CJTTING. 

VLLV(J)   =    PERCENT   OF    PREVIOUS  OLEVtll    TO   BE  LEFT    AT   REGN [ J) f 

INF ESTEO  STANOS. 
COMMON  /BLKA/BA,08H,DEN,HT,B0F,CFM,TCF 


COMMON  /B L KO /C BA S(2, 20), T BA S(2, 20), UBAS12, 20), C8F12, 201, T8F (2, 201 
1UBF(2,20I  ,COEN(2(20I ,T DEN (2,20  I , UDEN ( 2 , 2 0 t , CMCF ( 2 , 20 ) ,THCF(2,20»  , 
2MCF(2,20> , CTOT (2,20  1 ,TTOT (2,20  » , UTOT (2,20) , TAGE( 2,201 , U AGE ( 2,201  , 
3OM(2,20 I , UOM(2.20I  *  TH  T (2 , 20 1 , JHT (2,20) »  UDMR (2,20  1  , TOMR (2,201 

COMMON  /8LKE/A00HT(21 , AGE 0 ( 2 1 ,  S A SO ( 2 1 ,0ELAY(21 ,DLEV(2I,OMRT(2l, 
1HTSO (21 ,HTST (21 , 1 OPT (2) , JCYCL (21 , JO EN (2)  , KOL , K ST EP ( 2  I ,R EGN ( 3 1 , 
2RINT ,ROTA  *  TEMHI 2) , T IME  ,VLLV (2  I  , ROTAOS 

COMMON   /BLKF/IKC2) , KOUNT ( 2 ) , S3 AS , KPER , KO EL  1 2 ) 

NSTNO  =  0 
DO  10  J=16,18 
VAR(J)    =  0.0 
10  CONTINUE 


READ    (5,20)    NSTND, RINT 
20   FORMAT  (I5.F5.0) 

IF ( NSTND   .LE.   0)   GO  TO  120 
VARU6)    =  RINT 

C    READ   MINIMUM  COMMERCIAL  CUTS   FROM   CARD    TYPE  TWO. 

C 

RE AO    (5,30)   COMCU, COM6F 
30    FORMAT  (2F5.0) 
VARU7)    =  COMBF 
VAR(18)    =  COMCU 


C    EXECUTE    PROGRAM  ONCE    FOR   EACH  STAND. 


C    CHECK  FOR   UNWANTED  ZEROS    OR   BLANKS    IN  OATA. 
OO  <»0  1=1,7,2 

IFIVAR(I)    .LE.   0.0)   GO  TO  90 
t.0  CONTINUE 

00  50  1=13,19 

IF  ( V AR (  I)    .LE.   0.01    GO  TO  90 
50  CONTINUE 

C   COMPUTE   VALUES   FOR  FIRST   ENTRIES   IN   YIELD  TABLES. 


C   ENTER  LOOP   FOR  REMAINING  COMPUTATIONS  FOR 
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LIN  =    (ROTA  -  AGEO(l))    /    R I  NT  * 
IF IJDEN (2 )    . GT*    Ot    LI H  =  (?0TA 
□Q  70  K=1,LIM 
KPER  -  K 
CALL  CUTS 

S8AS   =   BAST(l)    ♦  BAST  1 2) 
CALL  PROJ 

SBAS   =   8AS0(1>    ♦  BASG(2) 
70  CONTINUE 

WRITE  YIELD  TABLES  FOR  A  STAND. 

00  80  1=1*2 

IF (A  GEO ( 1 1    .LE.    0.0)   GO  TO  80 
MOL    =  I 
CALL  TABLE 
80  CONTINUE 
GO  TO  110 


90    WRI T  £    (6,100)  STNO 
100   FORMAT    (1H1,///,10X,73HEXECUTION  STOPPED  BECAUSE   OF   NEGATIVE  OR  ZE 

IRQ  ITEM  ON  OATA  CAROS  OF  STANO,Fi».0) 
110  CONTINUE 

GO  TO  lt.0 
120    WRITE  (6,130) 

130   FORMAT    ( 1 H 1 , / // / , 3 0 x , 5 2H£ X£CU T I O N  STOPPED  BECAUSE   OF    ILLEGAL  VALUE 

1   OF  NSTND .  ) 
l*i0  STOP 
END 


UHTfl.J)    =  0.0 
UMCF(I,J)    =  0.0 
UTOTII.J)    =  0.0 
"*Q  CONTINUE 

OMLEV    =  0.0 
ICUT   =  0 
STNO   =  0.0 
TIME   =  0.0 

;    RE AO   TEST    DESCRIPTION  FROM  CARD   TYPE  THREE. 


C    RE AO    INTERMEDIATE  CUT    INSTRUCTIONS    F  R01    CARO    TYPE  F3UR. 

READ    (5,6  0  )    ICU  T, OMLE V .OELA Y( 1 ) , IOPT (1)  , THIN( 1 ) ,OSTY  (1) , JCYCL (1)  , 
10ELAY(2),  IOPT(2),TMIN(2)  ,  QSTY ( 21  ,JCYCL(2) 
60   FORMAT    ( 1 5  ,F5  . 1 .  2  ( F  5.  0  , 15 ,  2 F5 .  0 ,  15 )  ) 
ICUT  =    ICUT   *  1 

3   STORE   INITIAL    VALUES  FOR  ERROR  CHECKS  AND  LATER  USE. 
C 

VAR(7)    =  JCYCLU1 

IF(JCYCL(2)    .LE.    0)    JCYCL(2)    =  JCYCLd) 
VAR(8)    =  JCYCL(2» 
VARUM    =  THIN(l) 

IF(THIN(2)    .LE.    0.0)    THIN(2I    =  TMIN(l) 
VAR(15)   =  OSTY(l) 

IF(DSTY(2)    .LE.    0.0)    DSTr(2)    =  DSTY(l) 

Z    READ   INITIAL    STAND  VALUES   FROM  CARD    TYPE  FIVE. 

C 

RE AO    (5,70)  STnD,SIT£,AG£0(1»,D8H0(1>.HTS0(1>,QEN0(1),0MR(1I, 
1PINF (1 )  ,ST ART (1 ) , A GEO (2) , 03HO(2)  ,HTSO (2 ) ,DEN0(2)  , DMR  IZ) ,PINF(2)  , 
2START(2) 

70   FORMAT    ( 2 F 5. 0 , 2 ( F5 . 0 , 2F5 . 1 , F5 . 0 , 2F5  . 1 ,  F 5 .  0  )  ) 
;    STORE    INITIAL    VALUES    FOR   ERROR  CHECKS  AND   LATER  USE. 


SUBROUTINE  BEGIN 
TO    INITIALIZE    STAND    VARIABLES   AND    READ    IN    STAND  DATA. 


COMMON  /BLK0/C8AS(2 ,20 >,T9AS(2,20),USAS12»20>,CBF(2,20),TBF(2,20>, 
1UBF(2,20)  , COEN(2,20),TDEN(2.20>, UOENI2.2QI ,CMCF( 2,2  0)  ,TMCF(2,20) ,1 
2MCF(2,20)  ,CTOT(2,20),TTOT(2,201,UT0T(2,2Q) ,TAG£(2,20 ) ,U  AGE (2*20  I , T 
3DM(2,2Q),UOM(2,2Q),THT[2,2Q),UHT{2,2Q),UOMR(2,20>,TDMR(2,20> 

COMMON  /8LK£/ADDHT(  2)  ,  AGEO(  2)  *  BA  SO  (  2  )  ,  DEL  AY  (2  )  ,  )L  EV ( 2  >  *  DMRT  (2)  , 
1HTSO (2) ,HTST (2) ,I0PT(2) * JCYCL (21  * J0EN(2» . KOL, KS T EP ( 2 ) . REGN ( 3 ) , 
2RI NT, ROTA, TEMH(2),TI ME, VLLV(2>, ROTAOS 

COMMON  /BLKF/IK ( 2 1 ,K0UNT(2> , S 3 A S , KP E R , KOEL ( 2 ) 

DO  10  J=1.15 

VAR(J)    a  0.3 
10  CONTINUE 

DO  2  0  1=1,2 

AODHT ( I >    =  0.0 

AGEO(I)    =  0.0 

BASOU)    =  0.0 

BAST(I)    =  0.0 

08HO(I)    =  0.0 

DBHT(I)    =  0.0 

DELAY { I )    =  0.0 

OENO(I)    =  0.0 

OENT(I)    =  0.0 

DLEV(I)    =  0.0 

OMR(I)    =  0.0 

OMRT(I)    =  Q.O 

OSTY(I)    =  0.0 

HTSOdi    =  0.0 

HTST(I)    =  0.0 

IOPTIII    =  0 

IREGN(I)    =  0 

ISANS(I)    =  1000 

JCYCL(I)    =  0 

JDEN(I)    =  0 

KOUNT(I)    =  1000 

KSTEP(I)    =  1 

KTR(I>    =  0 

PINF(I)    =  0.0 

RESTU)    =  0.0 

S8D(I)    =  0.C 

SHC(I)    =  0.0 

START ( I J    =  0.0 

STF(I)    =  0.0 

THIN(I)    =  0.0 

VLLV(I)    =  0.0 
20  CONTINUE 

OO  30  1=1,3 

REGN(I)    =  0.0 
30  CONTINUE 

OO  <»0   1  =  1,2 

DO  i*0   J=l  .20 

CBASU.J)    =  0.0 

CBF(I,J)    =  0.0 

COENd.J)    =  0.0 

CMCF(I.J)    =  0.0 

CTOT (I.J)  =0.0 

TAGEfX.JI    =  0.0 

TBASU.J)    =  0.9 

T8F( If Jl    =  0.0 

TOEN(I.J)    =  0.0 

TOM(I,J»    =  0.0 

TOMR(I,J)    =  0.0 

THT(I.J)    =  0.0 

THCFCItJl    =  0.0 

TTOT (I, Jt    =  0.0 

UAGE(ItJ)    =  0.0 

UBASU.J)    =  0.0 

UBFtl.J)    =  0.0 

UQENUtJ)    =  0.0 

UOM(I.J)    =  0.0 

UOMR(I.J)    =  0.0 


VAR(l)  : 
VAR(2)  = 
VARI3I  = 
VAR(it)  = 
VAR(5)  - 
VAR16)  = 
VAR(9)  = 
VAR( 10 ) 
VAR( 11) 
VAR( 12) 
VAR( 13) 
VAR(21) 
VAR(22) 
JOEN(l) 
JOEN (2) 


OGEO(l) 
A  GEO  (2) 
OBHOC II 
DBHO  (2 ) 
OENO(l) 
D£NO(2) 
STARK  1) 
s  START ( 2  > 
:  OMR  (1 ) 
:  OMRI2) 
:  SITE 

:  PINF(l) 

■  PINF12) 
=  OENO(l) 
:  OEN012) 


READ  REGENERATION  CUT   INSTRUCTIONS   FROM    CARD   TYPE  SIX. 

RE AO    (5,80  )   REGN(l) ,  VLLV  ( 1  >  ,  REGN  ( 2  >  ,  VLL  V  12  >  ,REG-i  (  3) 
80   FORMAT    ( F 5  . 0 , F5 . 3 , F 5 . 0 , F5 . 3 , F5 . 0 ) 
VAR(19)    =  REGN(l) 
RETURN 
END 


SUBROUTINE  FIRST 

C 

Z    TO  COMPUTE    VALUES  FOR  FIRST   ENTRIES    IN    YIELD  TABLES. 
C 

COMMON    /BLKA/BA ,08H  ,0 £ N ,  H  T  ,  80 F  , C FM  ,  T C F 

C 

COMMON  /BLK8/8AST( 2) , DBHO (2) , OBHT ( 2 )  ,0EN0(2) , DENT (2)  .OMR  I  2)  ,PRET  , 
1REST (2) ,1 TM,M0L*DMLEV,IREGN(2) ,PINF (2) ,  ISANS<2> 

COMMON   /8LKC/0STY(2> ,KTR(2) ,ST NO, S30(2>, SITE, SMC (2) ,STF(2I  , 
1  THIN  ( 2 ) »  V  AR ( 22)  , C OMOF , CO MCU , S T AR T ( 2  I ,N0TE(16),i;UT 

COMMON  /BLKD/C8AS 12, 20 > ,TBAS( 2.20) ,U8AS 12.2  0) ,C3F (2,20)  ,TBF(2,20). 
1UBF(2,23) , COEN(2,2D) ,TDEN(2.20) , UOE N ( 2 . 20 ) • CMCF ( 2 , 2 0 ) . T MCF ( 2, 20  I ,U 
2MCF(2,20) ,CTOT(2,20),TTOT(2,20),UTOT(2,20),TAGE(2,20),UAGEI2,20),T 
3OH(2,20» ,UOM(2,20),THT(2,20).UHT(2,20),UOMR(2,20) ,TDMR(2,2Q) 

COMMON  /BL  KE/AODHT ( 2 ) , AG E O ( 2 ) , B A  SO ( 2 )  , QE L A Y ( 2 > , JL E V ( 2 > , OIR T ( 2 ) , 
1HTS0(2),HTST(2),I0PT(2),JCYCL(2I,J0EN(2),K0L,KSTEP(2),R£GN(3>, 
2RINT ,ROTA , T£MH(2) , T I ME , VLL V ( 2 ) , ROT AOS 

COMMON    /BLKF/IK (2» ,KOUNT(2),S8AS,KP£R,<0£L(2) 

KOL    =  1 
IRINT   =  RINT 

C    IDENTIFY    FINAL    AGE    IN    YIELD  TA3LE. 

OO  10  NA=i,3 
L  =  *»  -  NA 

IF (REGN (L 1    .GT.    0.0)   GO  TO  20 
10  CONTINUE 
20   ROTA  =   RE  GN( L  > 

VAR(20)    =  ROTA 

ROTAOS   =   AGEO(l)    ♦  OELAY(l) 

IF (ROT AOS    ,L T .   ROTA)   ROTAOS  =  ROTA 

Z    COMPUTE   VALUES   OF  STANO  CHARACTERISTICS,    ONE   STORY    AT    A  TIME, 


IF (OMR  II)    .GT.   Q.O)    GO   TO  70 
IF(PINFd)    ■  LE.    0.0)    GO  TO  <t0 
OMR! 1 1    =    10.0   »•    (-0.9^652   *   0  .  01651 
GO  TO   7  0 
t»0    TIME   =   AGEO(I>   -  ST  ART  ( 1 1 
I F  (  T  I  ME   .LE.    0.0)    GO   TO  100 
IF ( I    .EQ.    2)    GO    TO  60 
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□  MRU)    =    0. 120055   •    TIME   -   0.961031    •   T I H£  / 
1  TIME  •  SITE 
CO  TO  70 

60   [FtOENOIZl    > GT.    3000.0  )    GO  TO  65 

ONRI2I    =   0  .0  87923   •   T I  ME  -  0  .  0  000  18    •    TINE  • 
1  TINE   •  DNR(1>   *  OMR ( 1 )    •  OENOU) 
CO  TO  7  0 

65  OMR(2>    =    0.0<»M7S7  •   TINE  -  0.  00000  31.   •  TIH£ 
I  0.0  00051    •   TIME   •    OMR  (11    •  OHR(l)    •  OENO(l> 

70  IF(OMR(It  .L  T.  0.0)  OMR(I)  =  0.0 
IF  (OMR  ( I )  .GT.  6.0)  OMRU)  =  6.0 
IFIOHRUJ    .GT.    0.0)    ST  ART  ( I  t    =  1OO0Q.O 

C  COMPUTE  INITIAL  AVERAGE  HEIGHT   IF   NOT  REPORTED. 

100   IF(HTS0(I>    . GT.   0.0)   GO  TO  130 
IMAGEOd)    .GT.    55.0)    GO   TO  110 
MTSO(I)    =    Q.OH.i.1   •   AGEOII)    •  SITE 


0.1216?  •    AGEOd)   -   1.  50953 


GO  TO  120 
HO  HTSO(I)  = 
1  20.  52528 

HT50II)  =  10. 
120  MTPCT  =  1.073 

IFIHTPCT 

HTSOIX1 


5991.7  -  61.5019  /   AGEO(I)    »   0.80522   ■  ALOG10(SIT£» 
ALOGIO(SITE)    /  AGE 0 ( 1 ) 
.0  ••  NTSO(I) 

0.0367  •  OMR(I) 
0)    HTPCT   =  1.0 
HTSO(I)    •  HTPCT 


COMPUTE    INITIAL  8ASAL    AREA    A  NO  VOLUMES. 

130  84S0II)    =   OENO(I)    •   0.00Sit5<»2   •  OSHO(I) 
BA  =   B ASO  <  I) 
OBH  =   09H0  (I ) 
OEN  =  OENO(I) 
HT  =  HTSO(I) 
CALL  SMVOl 

*  STORE  VALUES   FOR  PRINTING. 

UAGEd.KOL)   =  AGEOII)    *  0.5 
MAX  =  OMR  I  I)    *    10.0    ♦    0. 5 
UONRII.KOL)    =  MAX 

UOHR(I.KOL)    =   UDMR (I »  KOL )    *  0.1 
UDEN(I.KOL)   =  OENO(I) 
MAX  a   BASO(I)    •  0.5 
UBASd.KOL)    =  MAX 
UOMII,KOL)    =  OBKO(I) 

MAX  =  -usum  ♦  0.5 

UHTfl.KOL)    =  MAX 
MAX  =    (TCF   •   0.1)    *■  0.5 
UTOT (I ,<OL I    =   MAX    •  10 
MAX  =    (CFM  •    0.  1)    *  0.5 
UHCF(I,KOL>    =   MAX    •  10 
MAX  =    C60F    •    0.01)    ♦  0.5 
UBF(I,KOL>    =   MAX    *  10C 

;  COMPUTE  VALUES  OF  NEEOEO  CONTROL  VARIABLES. 
C 

IK(I)    =   JCYCL(I)    /  IRINT 

ITEM  =    IRE&NI  1)    -   AGEOII))    /  R  I  NT   ♦  0.5 
IF  C I  X  C I >    .GT.    ITEM)    IK(I)    =  ITEM 
KD£L(I)    =   OEL  AY  ( 1 1    /  RINT   ♦  0.5 
REST(I)    =  THIN(I) 
OLEVdl    =  OSTTCI) 
lfcO  CONTINUE 
RETURN 
END 


IF (AGE 0 ( I )    . L T .   REGN (  1 ) )    GO  TO  20 
00  10  L=l 1 2 

IF (AGE 0 ( I >    .NE.   REGN(L))   GO  TO  10 
OLEV(I)    =    OLEV(I)    *  VLLV(L) 
REST(I)    =    OLE  V( I ) 

JCTCLd)    =  REGN ( L» 1 1    -  REGN(L)    ♦  0.5 
IK(I)    =  JCYCL(I)    /  IRINT 
IREGN(I)    =  1 
GO  TO  20 
10  CONTINUE 

COMPUTE  CURRENT  GROHING  STOCK  LEVEL  TO  DETERMINE  IF  STORY  IS  ALREADY 
AT   OR  BELOW  0 ESIRED   STOCKING  LEVEL. 


D3H0 ( I H 
QdHO  (  I) 


20  IF(D8H0(I(  . G 
GSLEVL  =  80.0 
GO  TO  CtO 

22    IF(09H0(I)    .GE.   5.0)  GO  TO  25 

GSLEVL    =   10.0  •   BASOUI    /    (11. 581.95 
GO  TO  it  a 

25    IF(08H0(I)    .GT.    10.0)    GO  TO  30 
TEH  =  OBHO(I)    •  DBMO(I) 

GSLEVL   =   80.0  •   BASO(I)   /    (7.76226  •   DBHO(I)    ♦   0.85289  »    TEH  - 
1  0.07952   •   TEM  •  DBHO(I)    -  3.'*562'*l 
GO  TO  <»0 
30    GSLEVL    =  9AS0d) 

;    CHANGE  O.6.H.    AND  O.M.R.    BY  PARTIAL  CUTTING. 

CO    IF  (OMR  ( I )    .LE.    0.0)    GO  TO  70 
IF ( I REGN ( I )    . EQ .    1)    GO  TO  70 
IFIDMR(I)    .GT.   DHL  E  V )    GO   TO  180 
IF(IOPT (I )    .GT.    0)    GO  TO  7] 
IF (KPER   .GT.    (KOEL(I)    ♦  1))    GO   TO  70 

C    COMPUTE    INITIAL   THINNING  LEVEL   IF    THINNING   FROM  A30VE.  EXECUTE 
Z    THINNING    FROM  A80VE. 

50    IF(THINII)    .LE.   60.0)    GO   TO  60 

REST(I)    =    THIN(I)    -    (OHR(I)    •    ITHIN(I)    -   60.0)    /  3.0) 

IF (OMR ( I )    .GT.   3.0)    REST  III    =  60.0 
60    IF (GSLE  VL    .LE.   RES  T  ( I )  )   GO  TO  180 

CALL  SHCUT2 

KSTEP(I)    =  2 

GO  TO  90 

C   EXECUTE   PARTIAL  CUTTING  FROM  BELOW. 

70    IF (GSL E VL    .LE.  REST(I))   GO   TO  180 
CALL  SHCUT1 

IF (PRE T   .GE.    100.0)    GO  TO  180 
KSTEP(I)    =  1 
90    IF ( B AST ( I )    . GE.    BASOdl)    GO  TO  180 

Z    COMPUTE   D.M.R.    AFTER  PARTIAL  CUTTING. 

C 

OMRT(I)    =    OMR ( I ) 

IF(OMRd)    .LE.    0.0)    GO  TO  110 

OELOMR   =   2.92759   *■   0.061CC  •  OMR(I)    *   OMR(I)    -  0.  2598ft  * 
1   SQRT(PRET)    -  0.00882   •  DMRII)    •   OMR(I)    •  SQRT I 3R£T ) 
IF(KSTEP(I>    .EQ.   1)    OELOMR  =   0ELDHR   •  0.66667 
IFtOELDIR    .LT.    0.0)    OELDMR  =  0.0 
ONRT(I)    =    OMR  II  I    -  OELOMR 
IFtOMRTd)    .LT.    0.0)    DMRT(I)    =  0.0 
IF(OHRT(I)    .GT.   6.0)    OMRT(I)    =  b.O 

C    COMPUTE   HEIGHT,    OAS AL    AREA,    DENSITY,    ANO   VOLUMES   AFTER   PARTIAL  CUTTING 


ALOG10 (PRET) 


SUBROUTINE  CUTS 

C  TO  CONTROL   EXECUTION  OF    INTERMEDIATE    ANO   REGENERATION  CUTS. 
COMMON    /8LKA/8A , D8 H , DE N, HT , 30 F , C FM, TCF 


COMMON  /8LKD/C8ASI2 ,20) , T9AS( 2,2  0) , US AS (2,2 0) , : 3F I  2 , 2 0 )  , T BF ( 2 , 2 0 ) , 
1U8F(2,2  0)  ,  CD EN (2, 20 ).T OEN (2, 20) ,UD£N(2,20) , CMCF ( 2 , 2 0  > , T MCF ( 2, 2 0 >  •  U 
2MCF(2,20)  , C TO T(2, 20), T TO TI2, 20), U TO T(2, 20), TAGE (2, 20), U AG E(2, 20), T 
30M(2,20>  ,U0M  (2,20)  ,THT  (2,20)  ,'JHT  (2,  20) ,UOMR (2,23 ) , T OMR (2. 20) 

COMMON  /8LKE/AD0HT (2) ■ AGEOI2)  ,  BA  SO ( 2) , J£LA Y ( 2 ) , 3 L E V ( 2 ) , DMRT I  2 1 , 
1HTS0(2) ,HTST 12) , I0PT(2) , JCYCL (2) , JO£N(2) , KOL , KSTEP ( 2 ) ,REGN(3) , 
2RlNT,ROTA,T£MH(2)  , T I M £  ,  V LL V  I  2 )  , ROTA  OS 


COMMON    /8LKF/IK(2> , KOUNT(2> , S 3 AS , KP £ R , KD£L ( 2 ) 


COMPUTE   VALUES  OF  STAND  CHARACTERISTICS,    ONE    STORY   AT    A  TIME. 
OO  230  1=1,2 

IF ( JDENt I )    .EQ.    0)    GO  TO  230 
IF ( I    .EQ.    1    .ANO.  AGEO(l) 
IF(I   . £Q.    2   .AND.  AGEO(2) 
BASTdl    =  BASOUI 
OBHT(I)    =  D9H0(I) 
DENT  d)    =    D£NO(  I ) 
OMRT(I)    =  DMR(I) 
HTST(I)    =    HT  SO ( I ) 
IF(KP£R   . LE.    KDEL1I))    GO   TO  230 


IF(KOUNT( I) 


K(I)  >   GO   TO  230 


"E0EFINE  INTERVAL  TO  FIRST  REGENERATION  ;UT  IF  LESS  THAN  CURRENT 
CUTTING  CYCLE. 

ITEM  =  (REGN(l) 
IF  (IREGH ( I >    . EQ. 
IREGN ( I I    =  0 
ITM  =  I 
XOUNT(I)    =  0 

TAGE  1 1 1  KOL  >   =  AGEOd)    *  0.5 
CHANGE  CUTTING  STANOARDS  IF  A   REGENERATION  CUT   IS  DUE. 


A  LUG 1 3 (PRET ) 


AlOGIO (BA3T (I ) > 


r  MCF ( I , KOL ) 


110    KS   =  KST£P(I) 

GO  TO  (120,130),  KS 
120    AODHT  (  I)    =    7. 61*833 

go  to  mo 

130   AOOHT(I)    =  2.C3837 

1   -  8.1.7960 
11.0    HTST(I)    =    HTSO(I)    ♦  AOOHT(I) 

JDENT    =    BAST  ( I )    /    (  0.  0  051.51.2    •    08HT  (  I )    *    D3HT(I)1    *  0.5 
DENT ( I )    =    JOE NT 

BAST(I)    =    0.  0051.51.2   *   D8HT(I)    •   OBHT(I)    •  DEN  T  d  ) 

BA  =   B  AST  (  I) 

OBH   =   □ 3HT  1 1  > 

OEN  =  OEN T  1 1 ) 

HT  =  HTST  (  I) 

CALL  SMVOL 

CHECK  TO   SEE    IF   INTERMEDIATE    CUT  SATISFIES   COMMERCIAL  CRITERIA, 
ANY.      DO    NOT    CUT    IF   COMMERCIAL  CRITERIA    NOT  SATISFIED. 
STORE   VALUES   FOR  PRINTING. 

MAX   =   CFM    •   0.1    ♦  0.5 
TMCF(I.KOL)    =  MAX  • 
IF ( TMCF (I , KOL)  .GT. 
CMCF (I, KOL  )    =  UMCF  t 
MAX    =   80F    •    0.01   ♦  0.5 
TBF ( I, KOL )    =   MAX   •  100 

IF (T  BF (  I, KOL)    .GT.    U6F  (I, KOL) )    UBF ( I , KOL )    =    T  8F ( I , KOL ) 
C8F(I,K0L)    =    UBF  < I , KOL )    -   T8F ( I , KOL ) 
IF( IREGN ( I  )    .fcQ.    1)    GO  TO  190 
GO  TO   (  17  0,  160,  190  ) ,  ICUT 
160    IF (KPER   . EQ.    (KO£L ( I )    »   1) >    GO  TO  190 
170    IF(C8F (  I, KOL)    .GE.    C0M8F)   GO  TO  190 

IF (CMCF ( I , KOL )    .GE.   COMCU)   GO   TO  190 
180    BAST ( I )    =  9ASO(I) 
09HT(I)    =  OBHOd) 
OENTd)    =  OENO(I) 
OMRT(I)    =  OMR(I) 
HTST(I)    =    HTSO ( I ) 
KTR(I)    =  1 
T  8AS ( I . KOL  ) 
TflFf I.KOLI 
TOEN ( I , KOL  ) 
TDM  t I, KOL ) 
TDMR ( I , KOL  ) 
THT ( I, KOL  ) 
TMCF ( I*  KOL) 
TTOT (I, KOL ) 
C3 AS (I.KOL) 
CBF( I, KOL ) 
COEN ( I , KOL ) 
CMCF ( I, KOL ) 
CTOT (I, KOL) 
GO  TO  230 
190    TOEN ( I , KOL )    =  OENT(I) 

MAX    =   OMR  T  ( I )    •    10.0    *  0 
TOMR ( I , KOL )    =  MAX 
TOMRd.KOL)    =   TOMR ( I , KOL 
MAX    =   HTST(I)    ♦  0.5 


U8AS( I, KOL) 
U8F  ( I , KOL ) 

UDEN ( I , KOL ) 
UO  M ( I , KOL ) 

UUMR( I, KOL) 
UHT (I.KOL) 

UMCF  (  I  ,KOL  ) 

UTOT  (I, KOL ) 


0.0 
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THT I  I , KOL  »  =  MAX 
MAX  =  BAST  (I)  ♦  0.5 
TBAS(I.KOL)  =  MAX 
MAX  =  TCF  •  0.1  ♦  0.5 
TTOT (  I, KOL)  =  MAX  •  10 
rOmi.KOL)  =  OBHT(I) 
COEN  I  I , KO  L  )  =  UOENII,KOL) 
ia.(I,«  .  i  =  UBASII.KOL) 
CTOTII.KOL)    =  UTOTII.KOL) 


TDEN ( I, KOL) 
T3 AS  1 1 . KOL I 
TTOT ( I , KOL ) 


IFIOMR(I)    .LE.   0.0)    CO  TO  195 

COMPUTE  PERCENTAGE  OF   TREES   INFECTEO  BEFORE   THE  PARTIAL   CUT ■ 

PlNFd)    =    57.  33  010    ♦  60.56935    *    AL O G  1  0  1 OMP I  I )  I 
IF  IP  INF [ I )    . L T .    0.0)    PINFd)    =  0.0 
IF(PINFd)    .  G  T.    100.0)    PINFtK    =  130.0 
IF(OMRTtl)    .LE.    O.C)    ISANSd)    =  0 

AOO   PERIODIC  CUTS  TO   TOTAL  YIELDS. 

195  STFII)    =   S  TF 1 1)    ♦  CTOT(I,<3L) 

IF (CMCF  1 1 , KOL)    .LT.   COMCU)   GO   TO  200 

SMCII)    =   SMCII)    ♦  CMCF(I,KOL) 
200   IF(C3F (  I, KOL)    .LT.    COMBF )   GO  TO  230 

S80<I>    =  SBDII)    *  C  8F  ( 1 1  <  OL  ) 

GO  TO  230 

REMOVE   STORY    I    IF  STORY    HAS  REACMEO   ROTATION  AGE. 

210  BAsrm  =  o.o 

09HUI)  =  0.0 

DENT ( I )  =  0.0 

OMRT ( I )  =  0.0 

HTST(I)  =  0.0 

8AS0(I>  =  0.0 

ODHO(I)  =  0.0 

OENO(I)  =  0.0 
DMR(I)    =  0.0 

HTSO(I)  =  0.0 

JOEN(I)  =  0 

ADD   FINAL   CUTS   TO   TOTAL  YIELDS. 

STF<I)    =    STF(I)    *  UTOTII.KOL) 

IFIUMCF  II , KOL)    .LT.    COMCU)    CO    TO  220 

SMC(I)  =  SMC(I)  ♦  UMCF(I,OLl 
220    IFIU8F 1 1 , KOL )    -LT.   COMBF)   GO  TO  230 

S3DII)  =  S3D(I)  ♦  U8F  1 1 ,  KOL  ) 
230  CONTINUE 

fit  TURN 

END 


GO  TO  i*0 

35   YOM  =   0 . 3  t«£»l  7  57"  -   0  .  00  00031.  •   OENT(2)    f   0.  000  051    •  DM 

1    OMR  T { 1 >    •  OENT(l) 
UQ    IFIOMRTII)    .GT ■    1.0    ■ AND •    DMRT(I)    .LT.    d . 5 )    GO   TO  50 

IF(I    •  EQ.    2  . ANO.   O ENT  (2 )    .GT.    3000  .  0)   GO  TO  5] 

OMR(I)    =    OHRT(I)    *■    0.5  »    YOM  *  RINT 

GO  TO  60 

50    OMR(I)    =    DHRT(X)    *    YOM   •  RINT 
60    IF(OMRil)    .LT.    OMR  T ( I ) )    OMR  ( I )    =  OIRT(I) 

IF(OMRII)    .LT.    0.0)    DMR(I)    =  0.0 

IF(OMPII)    .GT.    6.0)    OHR(I)    =  6.0 

IF(OMR(I)    .GT.   0.01    ST  ART  < I )    =    1000  0.0 

C 

;   COMPUTE  CURRENT  ACTUAL    GROWING  STOCK  LEVEL. 

70    I F ( D  AH T  ( I  ) 

GSLEVL   =  80.0 

GO  TO  90 
72    IFtOBHMIl  -GE 

GSLEVL    =  80.0 

GO  TO  90 

75   IF(OOMTII)    . GT.    10.0)    GO   TO  80 
TEH  =  DBHT(I)    *  DSHT(I) 

GSLEVL   =   80.0  *   SBAS  /    (7.7o226   •  OdHTIII    *■  0.45289 
12  *   TEM   ■    OBHT(I)    -  3.<*562i») 
GO  TO  90 
80   GSLEVL   =  SBAS 

Z  COMPUTE   NEW  O.3.H.   AT   ENO   OF  GROWTH  PROJECTION  PERIJO. 

90    IF (GSLEVL    .GT.    160.0)    GO  TJ  100 

09HOII)   =    1.0097  •   OBHT(I)    *■   0.009b  •   SITE   -  1.5766 
1   *  3.3021 
CO  TO  110 

100    08H0d»    =    0.88511    •    OBHTIII    ►    1.  29735    •    ALOGlu ( H  TST ( I >  >    *  0.00119 

1   •  08HMI)    •    SITE   ♦  62.37l7i»  /   S3AS   -  1.56975 
110    IF (DMRT ( I )    . LE*    3.5>   GO  TO  120 

OINC   =    (08H0(I)    -  03HT(D)    •    (  1.  197   -   0.056  •  J)MRT(I)) 

OflHOd)    =    DBHT(I)    *  OINC 
120    IOBHO  =   D 3H0 ( I )    *    10.0    *  0.5 

08HO (I )    =  IO0HO 

OBHO  (I)    =    OUHOl I )    •  0.1 

2   COMPUTE  CURRENT  NUMBtR   OF   TREES  »£R  ACRE   AND  BASAL  AREA. 


Tc M  -   J.  0  7  95 


ALOGi;(S1AS) 


)   GO  TO  II] 

0   .AND.  OM?T(II 

1*2271   •  DM?T(I1 


DIE   =  0.0 

IF  OMRT  (I I  .LT. 

IF(OENT(I)  .LT. 

OIE  =  20.  &6<»69  * 
1    ALO&10 (DENT ( I  I  I 

OIE  =  DTE   •  0.01 

IF(OIE    .LT.    0.0)    OIE   =  0.0 
130    OUT   =  0.0 

IF(D9HT(I)    .Ct.    10.0)    GO   TO  l'+O 

OUT  =  0  .  0  02*7  *  0.0012<t  •  JQHTU)  ♦  0. 
1   3.00000521  *   SdAS  *   SBAS  -  0.0000905 


IF(0UT 
1<*Q  IF(OIE 
JOENO  = 


.LT.  3.0)  OUT  =  0.3 
,  L  T .  OUT)  OIE  =  OUT 
OENT(l)    »    (1.0   -  OI 


DE  NO  (I )    =  JOENO 


SUBROUTINE  PROJ 
TO  MAKE   PERIOOIC  GROWTH  PROJECTIONS. 

COMMON    /3LKA/BA.03H  ,D£ N , HT , 30 F , CFH, T CF 


COMMON  /BLKO/C3ASI 2,2  0 )  ,TaAJ(2,20».U3AS(2.20).C3F(2.20).T9F(2,20), 
1UBF<2,23 1  ,  CD  EN  (2  ,20)  ,TOEN(2,20  )  ,  UOEN  (2,  20)  ,CMCF(  2,20  )  ,T'1CF  (2,23)  ,U 
2HCF(2,20) ,CTOT(2,20),TTOT(2,20),UTOT(2,23), T AGE ( 2,23) ,UAGE(2,23) ,T 
3OM(2.20),UOd2,20),THT(2,2J),JHT(2,20),JOMR(2,2J),TQMfc(2,20> 

COMMON  /3LKE/ AOUHT (2),AGE0(2),0AS0(2) , DELAY (2 ) , 
1HTSO  <2)  ,HTST (2)  ,I0PT(2)  . JCYCL ' 2)  ,JOEN  (2)  .KOL.KS 
2RINT,ROTA,TEMH(2),riME,VLLV(2).ROTAJS 

COMMON    /0LKF/IKC2)  ,  KOUNT  (2)  ,  S  3  A  S  ,  KP  ER  ,  KO  EL  (  2) 

KOL  =  KOL    ♦  1 

COMPUTE    VALUES   OF   ST ANO   CHARACTERISTICS,    ONt  STORY 

OO  220  1=1,2 

IF(JDEN(I)    . EQ.    0)    GO    TO  220 
AGEOd)    =    ACEO(I)    »  RINT 
KOUNT(I)    =    KOUNT II)    *  1 

COMPUTE    CURRENT    OWARF    MISTLETOE  RATING. 

ISANSd)    =    ISANSII)    •  1 
IFII SANSII >    .GT  .    2)    GO    TO  5 
FP  =  ISANSd) 
FP  =   0.025   *   FP   *  RINT 

OMR(I)    =    10.0   M    (-0.9U652   ♦   3.01651   *  (PINF(I) 
GO  TO  60 
5    IF(OHRT(I)    .GT.    0.0)    GO  TO  20 
OMR! I )    =  0.0 

TIME    =    AGEOII)    -    ST  ART  (I  I 
IFITIME    .LE.    0.0)    GO    TO  70 
IF  C I    .  E'J.    2   .ANO.  JDEN(l) 


1 .  50  95  3 


;    COMPUTE    CURRENT   AVERAGE  HEIGHT. 
C 

IF(GSL£VL    .GT ,    160.0)    GO   TO  200 
DO  190  J=l,2 
LUR   =  J 

GO  TO    1  15  0  ,  160) ,  LU3 
150    YARS    =  AGEO(I) 

GO  TO  170 
160    YArfS  =    AGE  OI I)   -  RINT 
170    IFIYARS   .GT.    55.0)    GO   TO  180 

T£MH(J)    =    0. 0 *    YARS    *    SITE    -    0.12162    *  YA*S 

GO  TO  190 

180    TEMH(J)    =    0.  599**7  -  61.5019  /   YARS   ♦   0.80522  •   A  LOG  10 ( S IT  El  ♦ 
1   20.52528    *  ALOG10ISITE)    /  TA*S 

TEMH1J)    =    10.0  TtMH(J) 
190  CONTINUE 

CHNG    =    TEMH(l)    -   T  tMH I  2  I 

GO  TO  210 

200   HTNEW   =   15.<*3021   ♦    1.107  •   HTSMI)    -   0.08637  •   AGEOII)    -   30-*. 12172 
1   /  SITE  -   0.  00021.1.7   •    SITE  •  SBAS 
CHNG   =    HT  NEW    -  HTST(I) 
210    PCT    =    l.C    -   0.  0002    •    OMRT  1 1 )    •    OMRT  I  I )    •    OMRT  d  I 
HTSO(I)    =    HTSTtI)    *■   CHNG    •  PCT 

C    COMPUTE    TOTAL    CU.    FT.    AND   CONVERT    TO   OTHER  UNITS. 

BA  =   9ASO  1 1) 
03H   =  08HOII) 
DEN  =  OENOII) 
HT   =  HTSOII) 
CALL  SWVOL 


STORE   VALUES   FOR  PRINTING. 

AGEO (  I) 


J  I    -.  .<    Ii  10 


U AGE (  I  ,  KOL ) 
MAX  =   OMR  ( I )    '   10.0    ♦  0.5 
UOMRII.KOL)    =  MAX 

UOMRII.KOL)    =   UOMRII.KOL)    •  0.1 
UOEN  d,  KOL)    =    DE  NO  (  I  ) 
MAX  =  BASO(I)    »  0.5 
U9ASd,K0L)    =  MAX 
UOM(I.KOL)    =  08MOII) 
MAX   =  HTSO(I)    ♦■  0.5 
UHT1I.OL)    =  MAX 
MAX  =   TCF    •   0.1    *  0.5 
MAX    *  10 


DMR(I)    =   0.120055   *    TIME   -   0.961031    •    TIME   /  OEHT(I)    -   0.00«*20  • 
1   TIME   •  SITE 
GO  TO  6  0 

10    IF(0ENT(2)    .GT.    3000.0)    GO    TO  15 

OMR(2)    =   0.087923   •   TIME  -  0.000018    •    TIME   •   OENTI2I    ♦   0.000051  • 
1   TIME  •  DMRT(l)    ♦   OMRT(l)    •   DENT  CI ) 
GO  TO  6  0 

15   OMRI2)    -   0.0<*<»1  757   •   TIME   -   O.OOQQ03<*   •   TIME   ♦   DENT  1 2 )  * 
1    0.  0  03051    •    TIME    •    DMRT I  1 )    •    OMRT(l)    •  OENT(l) 
GO  TO  60 

20    IFII    -EQ.    2   .ANO.   JOEN(l)    .GT.   0)    GO   TO   3 Cr 
IFIOENTII)    .LE.    15,0)    GO   TO  25 

YOH   =   0  .  1  200  55  -   0  .  961  031  /  OENTII)    -   0.  000  <*20   *  SITE 
GO  TO  <*0 
25    YOH    =  0.08 
GO  TO  <*Q 

30    IF(0ENT(2l    .GT.    3000.0)    GO   TO  35 

YOM  =  0.087923  -  0.  0O0018  •  OE NT  1 2 )  »  0.  000  05  1  •  UMRT(l)  •  OMRTI1) 
1   •   DENT  1 1 ) 


UTOTd.KOL) 
MAX   =   CFM  • 
UMCF II ,KOL 1 
MAX    =   90F  • 
UDFI I, KOL ) 
IF (KPEft  .GT 
220  CONTINUE 
RETURN 
ENO 


0.1 


0.5 
•  10 


0.01  *  0.5 
MAX  •  100 
KDELIII)    REST ( I ) 
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SUBROUTINE  TABLE 

TO  PRINT    YIELD    TABLES   A  NO    TABLE  FOOTNOTES. 

COMMON  /BLK3/BAST(2> ,OBHO(2), 03HT(2) , DENO(2) , OENT 12) ,QNR(2)  ,PR£T, 
1REST  121  1 1  T  Mt  HOL  ,DML  EV  ,  IR ECN  (2 >  ,  P  INF  <  2  »  ,  I  SANS  (  21 

COMMON  /BL  KC/DS T  Y  <  2  > »KTR 12) ,STnO,S80(2> • SITE.SM; (2)  ,STF(2) ■ 
1THIN (2 )  ,V ARI22)  ,C0M8F,C0MCU,START (2  I , NOTE ( 16) , ISUT 

COMMON  ✓8LK0/C8AS(2,2Q> , T9 AS  I  2 , 2 0  I , UB AS ( 2 , 2 3 > ,C3F(2,2Q),T9F(2,2QI 
1U8F1 2,23)  .  COEN ( 2,20 ) ,T0EN (2,20 1 ,UOEN (2,20) ,CMCF ( 2 , 2 0 ) , T MCF ( 2 ,2 0  I , 
2NCF(2.20),CTOT(2,20),TTOT(2,20),UTOT(2,2Q),TAG£(2,2Q)  ,UAG£ (2,2  0 )  , 
)OM(2 ,20) ,UOM (2,201  , THT 12,201  ,  JHT (2, 20) ,UOMR (2 ,20  > , TOMR ( 2* 20 ) 

COMMON  /SLKc/AOOHT (2) , AGfcO(2) , BASO( 2) , OtLAV C2 ) , 3LEV ( 2 ) * DNJRT 121 , 
1 hT  SO ( 2 ) ,HTST (2)  ,  IOPT(2  » , JCrCL(2) , JO EN  I  21 , KOl, KS TEP ( 2 ) ,R EGN ( 3 ) , 
2RINT,ROTA,TEMH(2),TlME,VLLV(2)  , ROT AOS 


280    FORMAT    (IH    , 10X,87HPRECOMH£RCIAL    INITIAL    THINNING    ALLOWED.  NONCOMM 
1ERCIAL    SUBSEQUENT    THINNINGS   NOT    ALL  OWED, ) 
GO  TO  310 
290   WRITE  (6,300) 

300    FORMAT    (1H    ,10X,35HALL    THINNINGS   MA  f    BE  NONCOMMERCIAL.) 
310    II    =    I    ♦  8 

IF(VAR(H)    .GT.    0.0    .AND.    VAR(II)    .  LT.    VAR(20))    GO   TO  330 
II   =    I    f  20 

IFlvAR(II)    .LE.    0.01   GO  TO  370 

WRITE    (6,320)    V AR ( I  I , V AR ( II I 
320    FORMAT    (1H    flQX,ldHAT    INITIAL    AG  E ,  F  5  .  0  ,  1 H  ,  ,  F5  .  0  ,  5  7H   PERCENT   OF  THE 
1    TREES    WERE    I NFECTEO    WITH  OWA3F  MISTLETOE.) 

GO  TO  37C 
330    WRITE    (6,31*0)    VAR ( 1 1  > 

3>*0    FORMAT    (1H    ,  10  X ,  <*0H  OWA  RF    MISTLETOE    INFECTION  STARTED   AT  A&E.F5.0) 
3  70    IF  (KTR  ( I )    .EQ.   0)    GO  TO  390 
WRITE    (6, 380 ) 

380    FORMAT    (1H    ,1QX,52HN0TE   THAT    NOT    ALL    SCHEDULED    THINNINGS    WE*E  PQSS 

LISLE.) 
390  RETURN 
END 


S  PRINT  HE  AO  I NG  S   FOR   YIELD  TABLE. 

If | I  .EQ.   2)   GO  TO  30 
IF(UOEN(2,l)    .GT.   0.0)    GO  TO  10 
WRITE    (6,5)  STNO 

5  FORMAT    (1H1,//, 32X, 60MYIELOS  PER  ACRE   OF   SOUTHWESTERN  PONOEROSA  PI 
1NE,  STAND   NUMBER, F7.0) 
GO  TO  20 
10  WRITE    (6,  15)  STNO 

15  FORMAT    (lHl,//t26X,60HYIELOS  PER  ACRE   OF   SOUTHWESTERN  PONOE^OSA  PI 

1NE,    STANO    NUMBER,F7.0,12H,    OVERS  TORT.  ) 
20   WRITE    (6,251    JS I TE , THI N ( 1 ) , OS T Y ( 1) 

25  FORMAT    I1HO,58X,10HSITE   INDEX  ,  I<*  /  H    ,  3  8  X  ,  29  H  T  H I  NNl  N&    INTENSITY--  I 
INITIAL-, F6.0,12H  SUBSEQUENT-, F6.0) 
GO  TO  5  0 
30  WRITE    (6,i«0)  STNO 

i.0  FORMAT  11H1,//,26X,  60HYIELOS  PER  ACRE  OF  SOUTHWESTERN  PONDEROSA  PI 
1NE ,    STAND    NUMBER, F7.0*13H,    UNOE  RSTORY . ) 

WRITE    (6,25)  JSITE,THIN(2),D5TY(2) 
50   WRITE   16,601    (NOTE  1  J ) , J= 1 , 16) 
60  FORMAT    1 1 H Q , 2 5X , 1 6 A 5 ) 

WRITE  (6,701 

70   FORMAT    (  1  HO  ,  2 1 X  ,  f.  1 H  CH  A  RA  C  T  E  RI  S  T  I  CS    BEFORE    ANO    4-TER    THINNING  ,  30  X  ,  2 
16HPERI00IC   INTERMEDIATE  CUTS) 
WRITE  (6,80) 

80   FORMAT  (1H0,<.X,5hSTANO,15X15h9ASAL.3X,7HAVERAGE,2X,7HAV£RAG£,3X,5H 
1 TOTAL  » 3X»9HMERCHANT-,3X, 9HSAWTIM8ER,9X,5HBASAL »tX,5HTOTAL»3X,9HMER 
2CHANT-, 3X.9HS AW  TIMBER) 
WRITE  (6,90) 

90   FORMAT      (1H   ,  5X,  3HAGE,  9X ,  5H  TREES  ,  3X  ,  <*H AREA .  <.  x  ,  6H  O  .  B .  H. ,  3  X  ,  6HHE I  GHT 
1,2X»6HV0LUME,2X,11HA8LE  V0LUM£ti,X,6HV0LUME,3x,5-(TKE£S,3X,t,HAR£A,3X 
2,6HV0LUME ,2X,11HA8L£    VOLUME, i.X»&HVOLUME) 
WRITE  (6,100) 

100  FORMAT  (1H  ,3*,7H(YEARS>  ,2X,3HONR,3X,3HNO.,3X,6HSQ.FT.,*.X,3HIN.,6X 
1.3HFT.  ,<«x  ,6HCU.FT.,5X,6HCU.FT.  ,  6X  .6H90.FT. ,  UX  ,  3  H  NO.  , 3X, 6HSQ .FT. , 2X 
2  .6HCU.FT. ,5X,faHCU.FT.,6X,6H8D.FT.) 

;  PRI NT  TABLE  ENTRIES  OF   STANO  CHARACTERISTICS. 

DO  150  N=l,20 

JAGEO  =  UAGE(I,N) 

IFtJA&EO   . EQ.   0)   GO  TO  120 

JOENO   =  UOEN(I,Nl 

JBASO  =  UBAS(I.N) 

OIAM   =  UOM(I.N) 

JHTSO   =  UHT(I,N) 

JTOTO  =  UTOT(I,N> 

JCFMO  =  UMCF(I,N) 

J90FO  =  U9F(I,N) 

WOMR    =  UDMR(I«N) 

WRITE    (6,110>    JAGEO, WOMR,J0EN0,JdAS0,DIAH,JHTS3, JTOTO, JCFMO, J DO FO 
110   FORMAT  (lH0,*.X,Ii,,2X,F5.1,2X,IS,2X,It.,5x,F5.1,5<,I3,'*X,I5,6x,I5,6X 
1,16) 

120    JAGEO   =    T  AGE ( I, N) 

IFl JAGEO   .EQ.   0)   GO   TO  150 

JOENT   =  TDEN(I.N) 

JBAST    =  T8AS(I,N) 

DNTR  =  TOH(I.N) 

JHTST    =  THT(I.N) 

JTOTT  =    T  TOT ( I, N) 

JCFMT    =  TMCF(I,N) 

JflOFT    =    T  BFI  I,N) 

JOENC  =  COEN(I,N) 

JBASC   =    CBAS  (  I,N) 

JTOTC   =  CTOT(ItN) 

JCFMC   =  CMCF(I,N) 

J9DFC   =  C3F(I,N) 

WDMRT  =  TOMR1I.N) 

WRITE    (6, 133)    JAGEO, WDMRT, JOENT, JBAST,UMTR, JHTST, JTOTT, JCFMT, J BO FT 
1  ,  JOENC, JBASC .JTOTC  »  JCFMC , J8DFC 
130  FORMAT    (1H   ,  U  X,  Id,  2  X,  F5.  1 ,  2X  ,  1 5,  2X ,  1 <•  *  5X  ,  F5  . 1 ,  5<  *  13.  **X.  15  , &X  ,  15  ,  6X 

1.16,  <»X,  15  ,  3X  ,  I3.5X,  m,  6X,  U,  8X(  15) 
150  CONTINUE 

JSTF  =  STF(I) 
JSMC  =  SMCII) 
JS80   =  SBO(I) 

WRITE    (6,160)    JSTF, JSMC, JSB3 
160    FORMAT    (1H0,/.67X,12HTOTAl    YI  £  L  O  S,  1  9X  .  I  5  ,  5X  ,  I  5  ,  ?  X  ,  1 6  ) 

e 

;  PRINT   APPROPRIATE   FOOTNOTES  FOR  TABLE. 
C 

WRITE  (6,170) 

170    FORMAT    (1H0,//, 1 1 X , 68 H MERCH .    CU.    FT.    -    TREES    6.0    INCHES    O.B.H.  ANO 
1  LARGER   TO   <♦  .  0-INCH  TOP.) 
WRITE  (6,180) 

180  FORMAT    (1H    ,10X,68HBO.   FT.    -   TREES    10.0    INCHES    D.B.H.   AND    LARGER  T 
10  VARIABLE    TOP   LIMIT. ) 
IFIIOPT(I)    .GT.   0)    GO   TO  200 
WRITE  (6,190) 

190  FORMAT    (IH    ■ 10X  ,67H  INITIAL   THINNING    FROM  ABOVE    ALLOWED    IN   STANDS  W 
1ITH  OWARF  MISTLETOE.) 
GO  TO  220 
200   WRITE  (6,210) 

210   FORMAT    (1H    , 10 X, 3 6H INITIAL    THINNING    MUST    8E    FROM    BELOW. ) 
220   WRITE    16  ,  230  )  OMLEV 

230   FORMAT    (IH    , 10X , 76HOWARF    MISTLETOE    RATING   ABOVE    WHICH  PERIODIC  THI 
INNINGS    WILL   NOT    8E    EXECUTED  -.Fd.l.lH.) 
WRITE    (6,2<tQ)  COMCU.COMBF 
2<.0   FORMAT    (IH    ,  10X,  <»<•  H  MINIMUM  CUTS    FOR    INCLUSION    IN    TOTAL    YIELDS--, F6 
1.0, 15H   CUBIC    FEET    ANO,F7.0,l2H    BOARD  FEET.) 
GOTO   (250,270,290),  ICUT 
250  WRITE  (6,260) 

260    FORMAT    I1H    ,10X,35HNO    NONCOMMERCIAL    THINNINGS  ALLOWED.) 

GO  TO  310 
270   WRITE    (6,  280  ) 


SUBROUTINE  SWVOL 

:    TO   COMPUTE    TOTAL   CU.    FT.    AND  CONVERT    TO   MERCH.    CU.    FT.    AND   BO ■  FT. 

COMMON  /9LKA/BA,0BH,0EN,HT,3DF,CFM,TCF 

FCTR  -  0.0 
PROO   =  0.0 

C    COMPUTE    TOTAL    CU .  FT. 


02H   =   03H    •  08H   *  HT 
IF (D2H    .GT.  5000.0)    GO  TO  10 
TCF   =   (  0.  53313  ♦  0  .  000  33   *  3A 
GO  TO  20 

10   TCF   =    (0. 00237  •   8A   ♦  0.00211 


0.00179  •  02H) 
D2H   -  1.093561 


DEN 
OEN 


.  0  -  INCH  TOP 


20    IF (DBH    .LT.   5.0)   GO   TO  70 

IF(08H    .GT.    6.5)    GO    TO  30 

FCTR  =   0.25222  •   DBH  -  1.01119 

GO  TO  50 
30    IF(08H    .GT.    1G.0)    GO   TO  HQ 

FCTR    =    3.Q2<*85  -    J. 09957   •   QSH   -    11.  3581*. 

GO  TO  50 

HQ   FCTR  =    1.  03936  -   l.t.103*.  /  03H 


VOLUMES    TO    VARIA9LE    TOP  IN 


50    IF(D8H    .LT.    8.0)    GO    TO  70 
IF  I03H    .GT  .    11.5)    GO   TO  60 

PROO  =  0.  0  028  •  BA  *  0. 0^355  •  OflH  *  DBH  -  2.78 
GO  TO  7  0 

60   PROD   =   0.8391.3   ♦   0.  20531   *  DBH 

CONVERT    TOTAL    CU.   FT.    TO  MERCH.    CU.    FT.    AND   90.  FT. 

70   BOF   =   TCF    »  PROO 
CFM   =    TCF    •  FCTR 
RETURN 
END 


SUBROUTINE  SWCUT1 


COMMON  /OL  KB/BAST (2)  ,D3HO(2) , 0  3HT (2) ,OENO (2),0E^T(2»,OMR(2) ,PR£T, 
1REST (2)  , ITM, MOL  »  DML  E  V ,  I RE GN ( 2 ) , P INF ( 2 ) , I S A NS ( 2 ) 


IF   OBHO    IS    LARGE  ENOUGH   FOR    BASAL    ARE  A    TO  REMAIN 


PRET    =  100.0 
OO  HQ  KJ=1,100 

IF (PRE  T   .LT.    50.0   .OR.    DBHO(I)    .LT.    3.0)   GO   TO  10 
08HE   =    3.73365  »   1.02008  •   08H0(I)    -   0.01107  •    (PRET   -  50.0) 
1  0.0001'*   •    (PRET   -  50.0)    •    (PR£T   -  50.01 
GO  TO  15 

10   POBHE   =  O.*t9*»01    *   0.71890   •  AL  OG 1 0  (  D8HO  (  I )  )    -    0.  2253  0  * 
1   ALOG1 0 (PRET )    ♦    0.12616  •   ALOS 1 3  1  OBHO  ( I ) I    *  ALJSIO(PRET) 
D8HE   =    10.0   »*  PD8HE 
15    IF ( I REGN ( I )   . EQ.   1)    GO  TO  20 
IF(OMR(I)    .LE.    0,0)    GO   TO  20 
TEM   =    D9HE   -   D8H0(I  ) 
OBHE  =   OBHO(I)   ♦  TEM  •  0.5 
OBHE     *   10.0   ♦  0.5 
IDBHE 

OBHE  •  0.1 

DENO(I)    •  PRET   *  0.01 

DENE   ♦  0.5 
NOENE 

0.  0  05t*5i*2  •    OBHE  •   38HE   •  DENE 

BASE   •   10.0    ♦  0.5 
NBASE 

BASE   *  0.1 

0.  0  05<*5<*2   •   OBHE  •  08h£ 


20  IOBHE 
OBHE 
OBHE 
OENE 
NOENE 
OENE 
BASE 
NBASE 
BASE 
BASE 
TMPY 


2? 


0.0)    GO   TO  190 

IF IT EM  .  LE . 
IF  t  T  EM  .LT. 

PRET   =   PRE T  -  1.0 

GO  TO  i*0 
30    PRET   =   PR E T  -  0.3 
40  CONTINUE 

GO  TO  140 

C 

C    COMPUTE    O.B.H.    IF  BASAL    AREA    INCREASES    WITH  D.3. 


50   PRET    =  40.0 

IFI08H01I1  .GT. 
00  1  30  J=  1. 100 
IF  C  PRET  . GE.  50 
PDBHE  =  0.49401 
1  ALOG10 ( PRET  >  * 
DBHE  =  10.0 
GO  TO  65 

60   OBHE   =   0. 73365  ♦ 
1   0.00014  ♦  (PRET 

65    IF (I REGN( I )  ,EQ. 
IFCOHR(I)    .LE.  0 


.AND.    03H0(I)    .GT.    2.9)    GO    TO  60 
0.71890  *   AL0G1Q(D3H0(I) >   -   0.22530  • 
0.12616  •  AL3S10 (DBHO (II »    *  ALOG 1C ( PRET 1 
PDBHE 

♦   1.02008  »  DBHQ  ( I)    -   0.01107  •    (PRET   -  5C 
-   50.01    ■    (P?ET   -  50.01 
11    GO  TO  70 
0)    GO  TO  70 


TEH 

OBHE  = 
70  IOBHE 
OBHE  = 
OBHE  = 
OENE  = 
NDENE 
OENE  = 
BASE  = 
NBAS  E 
BASE  = 
BASE  = 
BREAK 


DBHE 


.03682 


.10938 


OBHE  -  OBHO(I) 
=    03  HO ( 1 1    *   TEH  *  0 
=   DBHE    •    10.0    »  0. 
IOBHE 

OBHE  »  0.1 
OENO(I)    •    PRET  • 
=   OENE    ♦  0.5 
NOENE 

0. 0054542  *  OBHE 
=   BASE  *   10. 0    *  0.5 
N8ASE 

BASE  •  0.1 
=  49.9  *  REST  ( 1 1    /  8 
IF  ( BASE    .GT.   BREAK  1    GO  TO 
OBHP   =    (80.0   /  REST  (1)1  • 
GO  TO  100 
60   BUST   =   66.2  •  R£ST(I)    /  80.0 
IF  I  BASE   .GT.    BUST)    GO   TO  90 
DBHP   =    (80.0   /  REST (III 
GO  TO  100 
90    T MP Y   =   BASE  *  80.0   /  RES T ( 1 1 
TEH  =   TMPY  *  THPY 

DBHP  =  19.04740  •  THPY  -  0.26673  • 
1    -  448.76  833 

IF  ( T  MP  Y  .GT.  80.0)  OBHP  =  D8H0M) 
100    I08HP  =  OBHP  *   10.0    *  0.5 

OBHP    =  IOBHP 

OBHP    =    OBHP   •  0.1 

IFCOBHP  -   OB HE!  110.140,120 
110   PRET   =  PRET  •  1.02 

IF (PRE T   .GT.    100.01    GO  TO  190 

GO  TO  130 
120   PRET   =  PRET  *  0.98 
130  CONTINUE 
140    DBHMII    =  DBHE 

COMPUTE    POST-CUTTING  BASAL  AREA. 


IF(OBHMI)  .GT. 

SOFT   =  6.03633 

GO  TO  160 
145    IF(D8HT(I)    . GT. 

SQFT   =  11.58495 

GO  TO  160 
150    IF(OBHTd)  .GE. 

TEH  =   DBHT (I I  • 

SQFT   =  7.76226 
1   -  3.45624 
160   BAST ( I )    =    (REST (I) 

GO  TO  180 
170    BAST ( I  I    =    RES  T ( I  I 
180  RETURN 


10.0)    GO   TO  170 
DBHTII) 
DBHT(I)    ♦  0.85289 


T EM   -    0. 07952 


190    PRET   =  100. 
RETURN 
END 


SUBROUTINE  SHCUT2 


DUE    TO    THINHIN3  SOUTHWESTERN 


0  .OR.  OBHOdl  .LT, 
OSHOdl  •  0.  0  0807 
(PRET  -   50.0)    -  0. 


4.0)  GO  TO  10 
•  (PRET  -  50.0) 
91172 


I  =  ITM 

OETERMINE  PRET   AND  O.B.H.    AFTER  CJTTING   BY  AN  ITERATIVE  PROCESS. 

ITER  =  0 
PRET  P  =  100.0 
PRET   =  50.0 
5   ITER  =   ITER  ♦  1 
IF (PRET    .LT.  50 
03HE   =  0.98543 
1    (PRET   -  50.0) 
GO  TO  20 

10    PDBHE    =    0.51618    •    ALOG10(PRET(    t-    1.69219    •   ALOG1 0  (DBHO(  1 1  I 
1   -  0  .34768  •    ALOG1 0 (PRET 1    *   AL OG1 0 ( 03H0 1 1 ) I    -   1.  03421 

DBHE   =    10.0  •»  P08HE 
20    IF ( O HR ( 1 1    .LE.   0.01    GO  TO  30 

TEH  =  D3HE   -   08H0  < I  1 

OBHE   =   OBHOI I)    ♦  TEH  *  0.5 
30   OENE   =  OENOdl   •   PRET   •  0.01 

TREBA    =    0.0054542    •    OBHE    •  08HE 

BASE   =   TREBA   •  DENE 

IF  (DBHE  .GT.    2.0)    GO  TO  35 

GSLEVL    =   80.0   *   BASE   /    (6.03633  ♦  OBHE1 

GO  TO  60 
35   IF (OBHE 

GSLEVL 

GO  TO  60 

40    IF (OBHE   .GT.    10.076361   GO  TO  50 

TEH  =  OBHE  *  OBHE 

GSLEVL    =   80.0  •   BASE   /  (7.76226 
1  *  TEH  *   OBHE  -  3.4  56241 

GO  TO  60 
50   GSLEVL    =  BASE 
60   TEH  =  GSLEVL   -  REST(I) 

TEHA   =  ABS(TEM) 

IF (TEMA   .LT.    TREBA    .OR.    TEiA  .LT 
IF ( I TER   .GE.    100)    GO   TO  80 
OIFF   =   AB  S (PRETP  -  PRET)    •  0.5 
PRET  P   =  PRET 

IF ( TEM   .GT.   0.0)   GO  TO  70 
PRET   =   PRE  T   *  D IFF 
GO  TO  5 
70    PRET    =   PRET  -  OIFF 
GO  TO  5 


11.  097241 


0.011    GO   TO  3J 


SET    POST-CUTTING  VALUES. 

80    IDBHE  =  O0HE  *  10.0 
OBHE   =  IDBHE 
OlHTdl    =    D3HE  •  0.1 
IOENE  = 
DENT (I  I 
BAST  (I ) 
RETURN 
END 


OENE   *  0.5 

=  IOENE 

=  0.0054542 


22 


APPENDIX  2 

Output  of  Sample  Problem 


YIELOS  PER  ACRE   OF    SOUTHWESTERN   PONDEROSA   PINE,   STAND  NUMBER 


SITE   INDEX  70 
THINNING   INTENSITY —   INITIAL-     1Z0.  SUBSEQUENT- 
HYPOTHETICAL   HEALTHY   EVEN-AGED   STANO.    BEGIN  TABLE   AT   AGE  30. 
CHARACTERISTICS  BEFORE   AND  AFTER  THINNING 


PERIODIC   INTERMEDIATE  CUTS 


ST  AND 
AGE 
(YEARS) 

OMR 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
O.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BO. FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

30 
30 

0.0 
0 . 0 

950 
<>63 

119 
79 

<».a 

5.6 

25 
26 

1530 
940 

380 
380 

0 
0 

487 

40 

590 

0 

<>a 

0.0 

458 

109 

6.6 

35 

1520 

980 

0 

50 
50 

0.0 
0.0 

it  <.9 

226 

131. 

83 

7.4 
8.2 

44 
1.5 

2200 
1350 

1660 
1110 

500 
500 

223 

51 

850 

550 

60 

0.0 

225 

101) 

9.2 

51 

1870 

1640 

2200 

70 
70 

0.0 
0.0 

223 
1  39 

12* 
90 

10.1 
10.9 

58 
59 

2600 
1930 

2340 
1750 

5200 
5100 

84 

3", 

670 

590 

80 

0.0 

139 

107 

11.9 

64 

2560 

2360 

8400 

90 
90 

0.0 
0.0 

139 

89 

12<t 
90 

12.  8 
13.6 

69 
70 

3  22  0 
2350 

2990 
2200 

11200 
8500 

50 

34 

870 

790 

100 

0.0 

89 

103 

14.6 

74 

2890 

2730 

11100 

110 
110 

0.0 
0.0 

89 

28 

117 

45 

15.5 
17.2 

t3!f  ' 
80 

3440 
1360 

3260 
1310 

13800 
6000 

51 

72 

2080 

1950 

120 

0.0 

28 

53 

18.7 

83 

1680 

1620 

7900 

130 

0.0 

28 

62 

20.1 

86 

2010 

1950 

TOTAL 

10000 
YIELOS 

70  70 

5830 

MERCH. 
80.  FT. 
INITIAL 

CU.   FT.   -  TREES  6.0  INCHES  D.B.H 
-   TREES   10.0   INCHES   O.B.H.  AND 
THINNING  FROM   ABOVE   ALLOWED  IN 

.    AND   LARGER   TO   4.0-INCH  TOP. 
LARGER   TO  VARIABLE   TOP  LIMIT. 
STANDS   WITH   OWARF  MISTLETOE. 

OWARF   MISTLETOE  RATING   ABOVE   WHICH   PERIODIC   THINNINGS   HILL   NOT   BE   EXECUTEO   -  3.0. 
MINIMUM  CUTS   FOR   INCLUSION  IN  TOTAL    YIELDS--      320.    CUBIC  FEET   AND      1500.    BOARD  FEET. 
PRECOMMERCIAL   INITIAL   THINNING   ALLOHEO.    NONCOMMERCIAL   SUBSEQUENT   THINNINGS   NOT   ALL  OH  ED • 


SAWTIMBER 
VOLUME 
BO. FT. 


YIELDS   PER  ACRE   OF   SOUTHWESTERN   PONDEROSA   PINE,  STANO  NUMBER 


SITE   INDEX  70 
THINNING    INTENSITY--   INITIAL-     120.  SUBSEQUENT- 
HYPOTHETICAL   EVEN-AGED   STAND.    INFESTATION  STARTS   AT   AGE  10. 
CHARACTERISTICS   BEFORE   AND   AFTER  THINNING 


PERIODIC    INTERMEDIATE  CUTS 


STANO 

BASAL 

AVERAGE 

AVERAGE 

TOTAL 

MERCHANT- 

SAWTIMBER 

BASAL 

TOTAL 

MERCHANT- 

AGE 

TREES 

AREA 

O.B.H. 

HEIGHT 

VOLUME 

ABLE  VOLUME 

VOLUME 

TREES 

AREA 

VOLUME 

ABLE  VOLUME 

(YEARSI 

OMR 

NO. 

SQ.FT. 

IN. 

FT. 

CU.FT. 

CU.FT. 

BD.FT. 

NO. 

SQ.FT. 

CU.FT. 

CU.FT. 

30 

1. 

8 

950 

119 

4.  8 

25 

1530 

0 

0 

30 

6 

40  5 

41 

4.3 

23 

530 

0 

0 

545 

78 

1000 

0 

40 

1. 

0 

403 

74 

5.8 

32 

iooo 

450 

0 

50 

1. 

4 

399 

104 

6.9 

ill 

1620 

1120 

0 

50 

1. 

4 

3  99 

104 

6.  9 

41 

1620 

1120 

0 

0 

0 

0 

0 

60 

2. 

3 

353 

117 

7.8 

".7 

2020 

1600 

0 

70 

3. 

2 

299 

121 

8.6 

54 

2310 

1960 

1800 

70 

3. 

2 

2  99 

121 

8.6 

54 

2310 

1960 

1800 

0 

0 

0 

0 

80 

4. 

1 

242 

117 

9.<» 

60 

2500 

2200 

3500 

90 

4. 

9 

188 

107 

10.2 

64 

2500 

2260 

5100 

90 

4. 

9 

188 

107 

10.2 

64 

2500 

2260 

5100 

0 

0 

0 

0 

100 

5. 

4 

139 

92 

11.0 

69 

2310 

2100 

6300 

110 

5. 

8 

101 

78 

11.9 

72 

2080 

1920 

6800 

110 

5. 

0 

48 

45 

13.1 

73 

1220 

1140 

4300 

53 

33 

860 

780 

120 

5. 

4 

36 

41 

14.  5 

76 

1180 

1110 

4500 

130 

5. 

7 

27 

37 

15.9 

79 

1110 

1050 

",600 

TOTAL  YIELDS 

MERCH.  CU.  FT.  -  TREES  6.0  INCHES  D.B.H.  AND  LARGER  TO  4.0-INCH  TOP. 
80.  FT.  -  TREES  10.0  INCHES  D.B.H.  AND  LARGER  TO  VARIABLE  TOP  LIMIT. 
INITIAL  THINNING  FROM   ABOVE   ALLOHEO   IN  STANDS   WITH   OHARF  MISTLETOE. 

DHARF   MISTLETOE  RATING   ABOVE  HHICH   PERIODIC   THINNINGS  HILL  NOT    BE   EXECUTED   -  3.0. 
MINIMUM  CUTS  FOR  INCLUSION  IN  TOTAL    YIELDS--      320.    CUBIC  FEE  T   AND      1500.    BOARD  FEET. 
PRECOMMERCIAL   INITIAL    THINNING   ALLOHEO.    NONCOMMERCIAL   SUBSEQUENT   THINNINGS   NOT  ALLOHEO. 
OHARF   MISTLETOE   INFECTION  STARTED   AT   AGE  10. 
NOTE   THAT   NOT   ALL   SCHEDULEO   THINNINGS  WERE  POSSIBLE. 


SAWTIMBER 
VOLUME 
BD.FT. 
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YIELDS   PER  ACRE   OF    SOUTHWESTERN   PONOEROSA   PINE,    STAND  NUMBER 


SITE   INDEX  70 
THINNING   INTENSITY--   INITIAL-     120.  SUBSEQUENT- 

THO-STORI ED   INFESTED   STAND   WITH   IMHEOIATE   OVERSTORY   REHOVAL • 

CHARACTERISTICS  BEFORE   AND  AFTER  THINNING 


STAND 
AGE 
(YEARS) 


OMR 
<t.S 


TREES 
NO. 


BASAL 
AREA 
SQ.FT. 


AVERAGE 
D  .  3.  H. 
IN. 


AVERAGE       TOTAL       MERCHANT-  SAWTIMBcR 
HEIGHT     VOLUME      ABLE   VOLUME  VOLUME  TREES 

FT.  CJ.FT.  CU.FT.  BO. FT.  NO. 


PERIOOIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT-  SAWTIM8ER 
ABLE   VOLUME  VOLUME 
CU.FT.  BO. FT. 


TOTAL  YIELDS 


MERCH.  CU.  FT.  -  TREES  6.0  INCHES  O.B.H.  AND  LARGER  TO  <• . 0-INCH  TOP. 
BD.  FT.  -  TREES  10.0  INCHES  O.B.M.  AND  LARGER  TO  VARIABLE  TOP  LIMIT. 
INITIAL  THINNING  FROM   ABOVE   ALLOWED   IN  STANDS   WITH   DWARF  MISTLETOE. 

OWARF   MISTLETOE   RATING   ABOVE   WHICH   PERIOOIC   THINNINGS   WILL   NOT   BE   EXECUTEO   -  3.0. 
MINIMUM  CUTS  FOR  INCLUSION   IN   TOTAL   YIELDS—      320.    CUBIC  FEET    A  NO      1500.    BOARO  FEET. 
PRECOMMERCIAL   INITIAL    THINNING  ALLOWED.    NONCOMMERCIAL   SUBSEQUENT   THINNINGS   NOT  ALLOWEO. 


YIELOS   PER  ACRE   OF   SOUTHWESTERN   PONOEROSA  PINE,    STAND  NUMBER 

SITE   INOEX  73 
THINNING   INTENSITY--   INITIAL-     123.  SUBSEQUENT- 

THO-STORIED   INFESTEO   STAND   WITH    IMMEDIATE  OVERSTORY  REMOVAL. 

CHARACTERISTICS  BEFORE   AND  AFTER  THINNING 


3.,  UNDERSTORY. 


PERIOOIC   INTERMEDIATE  CJTS 


STANO 

BASAL 

AVERAGE 

AVERAGE 

TOTAL 

MERCHANT- 

SAWTIM8ER 

BASAL 

TOTAL 

MERCHANT- 

SAWTIMBEf 

AGE 

TREES 

AREA 

O.B. 

H. 

HE  IGHT 

VOLUME 

ABLE  VOLUME 

VOLUME 

TREES 

AREA 

VOLUME 

ABLE  VOLUME 

VOLUME 

YEARS) 

DMR 

NO. 

SQ.FT. 

IN. 

FT. 

CU.FT. 

CU.FT. 

BD.FT. 

NO. 

SQ.FT. 

CU.FT. 

CU.FT. 

BD.FT. 

30 

1 

0 

iOO 

55 

1). 

1 

21 

710 

0 

0 

30 

5 

512 

*Z 

3. 

9 

20 

560 

0 

0 

88 

13 

150 

0 

0 

1.0 

9 

509 

78 

5. 

3 

29 

1020 

330 

0 

50 

1 

I* 

503 

109 

6. 

3 

38 

1630 

91.0 

0 

50 

1 

it 

503 

109 

6. 

3 

38 

1630 

91*0 

0 

0 

0 

0 

0 

0 

60 

2 

3 

<t<t<t 

122 

7. 

1 

<."» 

2010 

1*5.0 

0 

70 

3 

2 

3  75 

123 

7. 

9 

51 

231.0 

1870 

0 

70 

3 

2 

375 

128 

7. 

9 

51 

231.0 

1870 

0 

0 

0 

0 

0 

0 

80 

«. 

P;  r 

30  3 

122 

8. 

6 

56 

2i>  <«0 

2060 

1900 

90 

k 

9 

231. 

110 

9. 

3 

61 

21.20 

2120 

3100 

90 

4 

9 

231. 

110 

9. 

3 

61 

21.20 

2120 

3100 

0 

0 

0 

0 

0 

100 

5 

3 

173 

96 

10. 

1 

65 

2280 

2050 

44  3  0 

110 

5 

8 

125 

81 

10. 

9 

69 

201.0 

1860 

5300 

110 

5 

0 

56 

1.5 

12. 

1 

70 

1150 

1070 

3800 

69 

36 

890 

790 

1500 

120 

5 

3 

1.3 

1.3 

13. 

5 

73 

1160 

1090 

1.200 

130 

5 

7 

32 

39 

11.. 

9 

76 

1100 

101.0 

1.300 

TOTAL 

YIELOS 

211,0 

1830 

5800 

MERCH. 

CU.  FT 

-  TREES  6. 0 

INCHES  O.B.H 

.   AND   LARGER   TO  4. 0-INCH  TOP. 

BD. 

FT 

.   -   TREES  10.0 

INCHES 

O. 

i.H.  ANO 

LARGER 

TO  VARIABLE  TOP 

LIMIT. 

INITIAL   THINNING   FROM   ABOVE    ALLOWEO   IN  STANDS   WITH   OWARF  MISTLETOE. 

DWARF   MISTLETOE  RATING   ABOVE   WHICH  PERIOOIC   THINNINGS   WILL   NOT   BE  EXECUTEO   -  3.0. 
MINIMUM  CJTS  FOR   INCLUSION  IN  TOTAL   YIELDS--      320.    CUBIC  FEET   AND      1500.    3JARD  FEET. 
PRECOMMERCIAL   INITIAL   THINNING   ALLOWEO.    NONCOMMERCIAL   SUBSEQUENT    THINNINGS   NOT  ALLOWEO. 
NOTE   THAT    NOT   ALL  SCHEOULEO   THINNINGS  WERE  POSSIBLE. 
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YIELDS   PER  ACRE   OF    SOUTHWESTERN   PONDEROSA   PINE,    STAND  NUMBER 

SITE    INDEX  70 
THINNING   INTENSITY--  INITIAL-     120.  SUBSEQUENT- 
TWO-  STORIED  INFESTEO   STAND.   OVERSTORY   REMOVAL   DELATED  20  YEARS. 
CHARACTERISTICS  BEFORE   AND  AFTER  THINNING 


PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
I  YEARS ) 

OMR 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CJ.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BO. FT. 

BASAL 
TREES  AREA 
NO.  SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

150 

<t.8 

20 

29 

16.1 

75 

800 

760 

3300 

160 

5.0 

16 

26 

17.1 

77 

71»0 

710 

3200 

170 

5.2 

12 

21 

18.0 

79 

630 

610 

2900 

SAWTIMBER 
VOLUME 
BO. FT. 


TOTAL  YIELOS 


MERCH.  CU.  FT.  -  TREES  6.0  INCHES  D.B.H.  ANO  LARGER  TO  4.0-INCH  TOP . 
BD.  FT.  -  TREES  10.0  INCHES  D.B.H.  AND  LARGER  TO  VARIABLE  TOP  LIMIT. 
INITIAL   THINNING  FROM   ABOVE   ALLOWED  IN  STANDS   WITH  OWARF  MISTLETOE. 

DWARF   HISTLETOE  RATING   ABOVE   WHICH   PERIODIC   THINNINGS  WILL   NOT   BE   EXECUTEO   -  3.0. 
MINIMUM  CUTS  FOR  INCLUSION  IN   TOTAL   YIELDS—      320.    CUBIC  FEET   AND     1500.    3OAR0  FEET. 
PRECOMMERCIAL   INITIAL   THINNING   ALLOWED.    NONCOMMERCIAL   SUBSEQUENT   THINNINGS   NOT  ALLOWED. 


STAND 
AGE 
(YEARS) 


YIELOS   PER  ACRE   OF    SOUTHWESTERN   PONDEROSA   PINE,    STANO  NUM3ER  l| 

SITE   INDEX  70 
THINNING   INTENSITY--  INITIAL-     120.    SUBSEQUENT-  90. 

TWO-STORIED   INFESTED   STAND.   OVERSTORY   REMOVAL   DELAYED  20  YEARS. 

CHARACTERISTICS  BEFORE   ANO  AFTER  THINNING 


UNDERSTORY . 


PERIODIC   INTERMEDIATE  CUTS 


BASAL        AVERAGE     AVERAGE       TOTAL       MERCHANT-        SAWTIMBER  BASAL         TOTAL        MERCHANT-  SAWTIMBER 

TREES        AREA  D.B.H.        HEIGHT     VOLUME      ABLE   VOLUME  VOLUME       TREES        AREA        VOLUME      A3LE   VOLUME  VOLUME 

NO.        SQ.FT.  IN.  FT.  CJ.FT.  CU.FT.  BO. FT.  NO.        SQ.FT.     CU.FT.  CU.FT.  BD.FT. 


30 

1.0 

600 

55 

4. 

1 

21 

710 

0 

0 

40 

1.5 

538 

70 

4. 

9 

29 

980 

0 

0 

50 
50 

2.5 
1.6 

471 
287 

83 
44 

5. 
5. 

7 
3 

38 
36 

1300 
680 

550 
220 

0 
0 

184 

39 

620 

330 

0 

60 

2.5 

253 

62 

6. 

7 

<t3 

1010 

670 

0 

70 
70 

3. <t 
3.4 

213 
213 

73 
73 

7. 
7. 

9 

9 

49 
49 

1290 
1290 

1030 
1030 

0 
0 

a 

0 

0 

0 

0 

80 

4.2 

172 

75 

9. 

0 

55 

1470 

1270 

1400 

90 
90 

5.1 
5.1 

133 
133 

74 
74 

10. 
10. 

1 
1 

60 
60 

1590 
1590 

1430 
1430 

3000 
3000 

0 

0 

0 

0 

0 

100 

5.5 

98 

66 

11. 

1 

64 

1530 

1400 

4200 

110 
110 

5.9 
5.7 

71 
50 

58 
45 

12. 
12. 

2 

9 

67 
68 

1430 
1140 

1320 
1060 

4800 
4000 

21 

13 

290 

260 

800 

120 

6.0 

36 

40 

14. 

2 

71 

1050 

99  0 

3900 

130 

6.0 

26 

35 

15. 

5 

73 

95  0 

910 

TOTAL 

3900 
YIELOS 

1860 

1240 

3900 

MERCH. 
BO.  FT 

CU.  FT. 
-  TREES 

-  TREES  6.0 
10.0  INCHES 

INCHES 
D.B.H 

D.B. 

AND 

H.  ANO 
LARGER 

LARGER   TO  4.0-INCH  TOP. 
TO  VARIABLE   TOP  LIMIT. 

INITIAL  THINNING  FROM   ABOVE   ALLOWEO   IN  STANDS   WITH   DWARF  MISTLETOE. 

DWARF   MISTLETOE  RATING  ABOVE   WHICH  PERIODIC   THINNINGS   WILL  NOT   BE  EXECUTED   -  3.0. 
MINIMUM  CJTS   FOR  INCLUSION  IN  TOTAL   YIELDS--      320.    CUBIC   FEET   ANO      1500.    80AR0  FEET. 
PRECOMMERCIAL   INITIAL   THINNING   ALLOWED.    NONCOMMERCIAL   SUBSEQUENT   THINNINGS   NOT  ALLOWED. 
NOTE   THAT    NOT   ALL  SCHEDULED   THINNINGS  WERE  POSSIBLE. 
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Agriculture— CSU,  Fort  Collins 
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